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ABSTRACT 



An experimental study of laminar flow heat transfer of an ethylene glycol/water 
mixture in an electrically heated horizontal tube using wire mesh (HEATEX) and 
twisted tape inserts was investigated. Twelve thermocouples, inserted in the tube wall 
at four longitudinal locations, enabled a mean inside experimental heat-transfer 
coefficient to be accurately measured. A constant wall heat flux boundary condition 
was placed on the wall by wrapping six 200 W flexible heater tapes tightly around the 
tube. The ethylene glycol/water mixture provided a coolant Reynolds number 
between 200-5000 and a Prandtl number between 30-140. Two smooth inside 
diameters and a roped tube profile were tested with and without the inserts. Heat- 
transfer correlations for tubes without inserts were developed and compared with 
theory for both thermally and hydrodynamically developing flow. Correlations were 
also developed for the two types of inserts. Nusselt numbers for fully developed flow 
were found to be a function of Reynolds and Prandtl numbers for the wire mesh 
insert and a function of tape twist ratio, Reynolds and Prandtl numbers for the 
twisted tape insert. Heat transfer enhancements of over 7 for the wire mesh insert 
and over 4 for the twisted tape insert at high Reynolds numbers were obtained over 
the empty tube. By using these correlations in conjunction with earlier obtained 
refrigerant condensation data (using the same tubes, inserts, and coolant), more 
accurate values of the outside condensation heat-transfer coefficients were obtained. 



in 



/ f f; V'- 

(jHc Slip 

o t 



TABLE OF CONTENTS 

I. INTRODUCTION 1 

A. BACKGROUND 1 

B. REFRIGERANT CONDENSATION RESEARCH AT NPS 2 

C. OBJECTIVES 7 

n. LITERATURE SURVEY 8 

A. INTRODUCTION 8 

B. LAMINAR FLOW IN CIRCULAR DUCTS 9 

1. Boundary Layers and Entrance Lengths 9 

a. Hydrodynamic Boundary Layer 9 

b. Thermal Boundary Layer 10 

2. Fully Developed Flow 12 

3. Entry Length Problem 14 

a. Thermal Entry Region - Uniform Wall Temperature ... 15 

b. Thermal Entry Region - Constant Wall Heat Flux 17 

c. Combined Entry Length 17 



IV 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOl 
MONTEREY CA 93943-5101 



C. AUGMENTATION TECHNIQUES TO ENHANCE TUBESIDE 

HEAT TRANSFER 19 

1. General Introduction 19 

2. Heat Transfer Enhancement with Twisted Tape Elements ... 20 

3. Heat Transfer Enhancement with Wire Mesh (HEATEX) 

Elements 22 

III. EXPERIMENTAL APPARATUS 25 

A. SYSTEM OVERVIEW 25 

B. TUBES AND INSERTS TESTED 28 

1. Tubes Tested 28 

2. Insert Elements Tested 29 

C. COOLANT SYSTEM 31 

1. Primary Coolant System 31 

2. Secondary Coolant System 31 

D. INSTRUMENTATION AND TEST SECTION 32 

1. Instrumentation 32 

2. Thermocouple Attachment Method 34 

3. Test Section 36 

E. DATA ACQUISITION AND CONTROL SYSTEM 37 

IV. EXPERIMENTAL PROCEDURES 39 

A. CALIBRATION EXPERIMENTS 39 



B. PROCEDURES FOR DETERMINATION OF INSIDE HEAT 

TRANSFER RATE 40 

1. General Procedures 40 

2. Specific Experiments 42 

a. Copper Smooth Tube 42 

b. Cu/Ni Externally Finned Tube 43 

c. Cu/Ni Korodense Tube . 44 

V. DATA REDUCTION 45 

A. AVERAGE INSIDE HEAT TRANSFER COEFFICIENT 45 

B. ADDITIONAL TUBE LENGTHS 53 

C. LOCAL NUSSELT NUMBERS 55 

VI. RESULTS AND DISCUSSION 56 

A. GENERAL 56 

B. SMOOTH TUBE 58 

1. Smooth Tube - No Insert 58 

2. Smooth Tube - Twisted Tape 61 

3. Smooth Tube - HEATEX Insert 61 

C. CU/NI EXTERNALLY FINNED TUBE 72 

1. Cu/Ni Finned Tube - No Insert 72 

2. Cu/Ni Finned Tube - Twisted Tape Insert 73 



vi 



3. Cu/Ni Finned Tube - HEATEX Insert 73 

D. KORODENSE TUBE 84 

1. Korodense Tube - No Insert 84 

2. Korodense Tube - Twisted Tape Insert 85 

3. Korodense Tube - HEATEX Insert 86 

E. INSERT COMPARISON BETWEEN TUBES 96 

F. LAMINAR FLOW ENHANCEMENT RATIO 101 

G. REPROCESSED REFRIGERANT CONDENSATION DATA 107 

vn. conclusions no 

A. GENERAL 110 

to 

B. SPECIFIC no 

VHL RECOMMENDATIONS 112 

APPENDIX A. FLOW METER CALIBRATION 113 

APPENDIX B. HEAT LOSS TO AMBIENT EXPERIMENT 116 

APPENDIX C. FRICTIONAL/HEAT LEAKAGE EXPERIMENT 120 

APPENDIX D. POWER CALIBRATION 122 

APPENDIX E. TUBE AND INSERT SPECIFICATIONS 125 

APPENDIX F. UNCERTAINTY ANALYSIS FOR THE INSIDE HEAT 

TRANSFER COEFFICIENT 127 

APPENDIX G. PROGRAM LISTING 141 



vn 



LIST OF REFERENCES . 156 

INITIAL DISTRIBUTION LIST 159 



viii 



LIST OF TABLES 



TABLE 3.1 DATA ACQUISITION UNIT CHANNEL ASSIGNMENTS . . 38 

TABLE 4.1 DATA FILE NAMES 41 

TABLE 6.1 SUMMARY OF CORRELATIONS 57 

TABLE 6.2 SUMMARY OF LAMINAR FLOW ENHANCEMENT 

RATIOS 102 

TABLE A.l. MASS FLOW RATE (KG/S) AT -15, 0 AND 24°C 114 

TABLE A-2. FLOWMETER CALIBRATION REGRESSION RESULTS 114 

TABLE B.l EQUILIBRIUM AVERAGE WALL TEMPERATURE 

RESULTS 117 

TABLE C.l HEAT LEAKAGE RESULTS 121 

TABLE D.l POWER CALIBRATION RESULTS 124 

TABLE E.1 TUBE CHARACTERISTICS/SPECIFICATIONS 125 

TABLE E.2 INSERT CHARACTERISTICS/SPECIFICATIONS 126 

TABLE F.l SUMMARY OF RESULTS FOR SAMPLE 

UNCERTAINTY ANALYSIS 140 



IX 



LIST OF FIGURES 



Figure 1.1 NPS Evaporator/Condenser Refrigeration Apparatus 4 

Figure 2. 1 Laminar Hydrodynamic Boundary Layer Development in a Circular 

Tube (Courtesy [Ref. 10]) 10 

Figure 2.2 Thermal Boundary Layer Development in a Heated Circular Tube 

(Courtesy of [Ref. 10]) 11 

Figure 2.3 Local Nusselt Number Obtained from Entry Length Solutions for 

Laminar Flow in a Circular Tube (Courtesy [Ref. 10]) 15 

Figure 3.1 Primary and Secondary Coolant Systems 26 

Figure 3.2 Data Acquisition and Control Systems 27 

Figure 3.3 Wolverine Korodense Tube (Courtesy [Ref. 21]) 29 

Figure 3.4 Twisted Tape and HEATEX Inserts 30 

Figure 3.5 Instrumented Tube 33 

Figure 4.1 Heat Balance 40 

Figure 5.1 Temperature Rise Corrections 48 

Figure 5.2 Cu/Ni Finned Tube Temperature Rise Corrections 51 

Figure 5.3 Additional Tube Lengths 54 

Figure 6.1 Smooth tube - No insert (Nu m vs. Re) results. 63 

Figure 6.2 Smooth tube - No insert (Nu m vs. Re) laminar flow results 64 

Figure 6.3 Smooth tube - No insert (Nu m vs. Re) transitional flow results. . . 65 

Figure 6.4 Smooth tube - No insert (Nu m vs. x*) results 66 



x 



Figure 6.5 Smooth tube - No insert laminar flow correlation results 67 

Figure 6.6 Smooth tube - Twisted tape (Nu m vs. Re) results 68 

Figure 6.7 Smooth tube - Twisted tape correlation results 69 

Figure 6.8 Smooth tube - HEATEX (Nu m vs. Re) results 70 

Figure 6.9 Smooth tube - HEATEX correlation results 71 

Figure 6.10 Cu/Ni externally finned tube - No insert (Nu m vs. Re) results. . . 75 

Figure 6.11 Cu/Ni externally finned tube - No insert (Nu m vs. Re) laminar flow 

results 76 

Figure 6.12 Cu/Ni externally finned tube - No insert (Nu m vs. Re) transitional 

flow results 77 

Figure 6.13 Cu/Ni externally finned tube - No insert (Nu m vs. x*) results. ... 78 

Figure 6.14 Cu/Ni externally finned tube - No insert laminar flow correlation 

results .— ..... 79 

Figure 6.15 Cu/Ni externally finned tube - Twisted tape (Nu,,, vs. Re) 

results 80 

Figure 6.16 Cu/Ni externally finned tube - Twisted tape correlation results. . 81 

Figure 6.17 Cu/Ni externally finned tube - HEATEX (Nu m vs. Re) results. . 82 

Figure 6.18 Cu/Ni externally finned tube - HEATEX correlation results. ... 83 

Figure 6.19 Korodense tube - No insert (Nu m vs. Re) results 87 

Figure 6.20 Korodense tube - No insert (Nu m vs. Re) laminar flow results. . 88 

Figure 6.21 Korodense tube - No insert (Nu m vs. Re) transitional flow 

results 89 

Figure 6.22 Korodense tube - No insert (Nu m vs. x*) results 90 

Figure 6.23 Korodense tube - No insert laminar flow correlation results. ... 91 



xi 



Figure 6.24 Korodense tube - Twisted tape (Nu m vs. Re) results 92 

Figure 6.25 Korodense tube - Twisted tape correlation results 93 

Figure 6.26 Korodense tube - HEATEX (Nu,,, vs. Re) results 94 

Figure 6.27 Korodense tube • HEATEX correlation results 95 

Figure 6.28 No insert comparison under laminar flow for all three tubes. ... 97 

Figure 6.29 Twisted tape insert comparison under laminar flow for all three 
tubes 99 

Figure 6.30 HEATEX insert comparisons under laminar flow for all three 

tubes 100 

Figure 6.31 Insert condition comparison for laminar flow in a smooth tube. 104 

Figure 6.32 Insert condition comparison for laminar flow in a Cu/Ni externally 
• finned tube 105 

Figure 6.33 Insert condition comparison for laminar flow in a Korodense 

tube 106 

Figure 6.34 Reprocessed refrigerant condensation data using new 

correlations 109 

Figure A.1 Flowmeter Calibration 115 

Figure B.l Heat loss to ambient; Average wall temperature versus time for 

1W and 7W ' 118 

Figure B.2 Heat loss to ambient; Average wall temperature versus heat • 
input 119 

Figure D.l Electrical schematic of power calibration experiment . 124 



xu 



NOMENCLATURE 



A 

D 

Gz 

h 

H 

I 

k 

L 

LMTD 

m 

Nu 

P 

Pr 

<? 

q" 

r 

R 

R w 



area, m 2 

specific heat, J/kg K 
diameter, m 
Graetz number 

heat transfer coefficient, W/m 2 K 
pitch for 180° rotation of tape 
current, amps 

thermal conductivity, W/m K 
length, m 

log mean temperature difference 
mass flow rate, kg/s 
Nusselt number 
perimeter, m or Power, W 
Prandtl number 
heat transfer rate, W 
heat flux, W/m 2 
radius, m 

electrical resistance, ohms 



xm 



wall resistance, K/W 



Re Reynolds number 

Re, Reynolds number for swirl flow 

T Temperature, K 

t fin thickness, m 

U overall heat transfer coefficient, W/m 2 K 

u velocity component in x direction, m/s 

v velocity component in y direction, m/s 

V velocity, m/s or Voltage, volts 

x’ dimensionless axial coordinate, (x/D)/(Re Pr) 

y ‘ twist ratio, H/D 

Greek Symbols 

8 tape thickness, m or boundary layer thickness 

€ heat transfer enhancement ratio 

li dynamic viscosity, kg/m s 

v kinematic viscosity, m 2 /s 

p density, kg/m 3 

Subscripts 

c coolant or cross section 

D tube diameter 

fd fully developed 

f ' fin 

h hydrodynamically 
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HL 


heated length 


i 


inside or inlet 


lam 


laminar 


m 


mean 
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outside or outlet 


s 


surface 


sat 


saturation 


t 


thermally 


t/c 


thermocouple 


X 


local 
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wall 
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I. INTRODUCTION 



A- BACKGROUND 

With the increasing complexity of shipboard weapons and combat systems used 
aboard United States Naval vessels, the need for more energy efficient, lightweight, 
and high capacity refrigeration, air conditioning, and chilled water systems has grown 
tremendously. One way to meet this need is to redesign the refrigeration condenser. 
In redesigning these condensers, changes in condenser tube material can be 
combined with various heat transfer enhancement methods to achieve a lighter 
weight, more efficient condenser tube. 

Emphasis on refrigerant condensation research reached a much higher level of 
importance due to the pressing environmental issues of ozone depletion and global 
warming. In September 1987, an international conference was held in Montreal, 
Canada, in which 24 nations representing the United Nation’s Environment Program 
(UNEP), signed the Montreal Protocol on Substances that Deplete the Ozone Layer. 
This agreement called for a near term freeze on the production and consumption of 
chlorofluorocarbons (CFCs), the major contributor to ozone depletion 
[Ref. 1]. In August 1988, the Environmental Protection Agency (EPA) 
adopted the regulations of the Montreal Protocol under the Clean Air Act [Ref. 2]. 
These regulations required a freeze in both production and consumption of CFC-11, 
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12, 113, 114 and 115 at 1986 levels. The U.S. Navy currently uses CFC-12 and CFC* 
114 in its shipboard air conditioning and refrigeration systems. In June 1990 at a 
progress meeting held in London, England, UNEP delegates agreed to even larger 
reductions, calling for a 50% reduction by 1995, 85% reduction by 1997, and for 
CFCs to be totally phased out by the year 2000 [Ref. 3]. In the spring of 1992, 
President Bush moved up the CFC phase out deadline to 1995. 

Because of the Navy’s commitment to the CFC phase out policy and the 
increasing need for lighter weight, more efficient refrigeration systems, refrigerant 
condensation research has surfaced to the forefront of Naval Research. 

B. REFRIGERANT CONDENSATION RESEARCH AT NPS 

In order to examine heat transfer characteristics of various condenser tube 
designs and alternative refrigerants, a condenser/evaporator test platform was 
constructed at the Naval Postgraduate School. Construction was begun by Zebrowski 
[Ref. 4] and completed by Mabrey [Ref. 5] in 1988. The test platform, shown in 
Figure 1.1, consisted of an evaporator and a condenser section with an associated 
ethylene glycol/water mixture coolant system. 

The evaporator section was designed so that in addition to supplying the 
refrigerant vapor to the condenser, it could also be utilized to conduct bundle boiling 
experiments. The condenser section consisted mainly of four instrumented condenser 
tubes and a series of auxiliary condenser coils used to control pressure in the 
apparatus. Details of design, construction and operation of the condenser/evaporator 
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test platform are provided by Zebrowski [Ref. 4], Mabrey [Ref. 5] and Mazzone 
[Ref. 6]. 

Unlike a typical vapor compression refrigeration system used aboard ships, the 
test platform at NPS lacked a compressor between the evaporator and the condenser. 
This compressor is normally required to increase the pressure (and hence 
temperature) of the refrigerant vapor sufficiently enough so that it can give up its 
latent heat to the ambient sea water, the primary coolant. Because there was no 
compressor in the NPS test platform, the condenser must operate at the saturation 
conditions of the evaporator. Therefore, a much colder coolant, a refrigerated 
ethylene glycol/water mixture, was needed to condense the vapor, leading to laminar 
flow conditions within the condenser tubes. This in turn leads to poor heat transfer 
due to a large inside thermal resistance. 

The primary objective of the condensation experiments of Mabrey [Ref. 5] and 
Mazzone [Ref. 6] was to accurately determine outside heat transfer coefficients. 
Briefly, this was accomplished by first determining the heat transfer rate to the 
coolant by: 

q = m c p A T (1-1) 

where m is the coolant mass flow rate and AT the coolant temperature rise. 
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THERMOCOUPLE ■*» 1 THERMOCOUPLE ** 2 




Figure 1.1 NPS Evaporator/Condenser Refrigeration Apparatus. 



The overall heat transfer coefficient, U 0 , could then be determined by: 



U = 

° A. LMTD 



( 1 . 2 ) 



where A„ is the outside surface area of the tube and LMTD is the log mean 
temperature difference, defined as. 



LMTD = 



The overall heat transfer coefficient can also be thought of in terms of a sum of 
thermal resistances. The overall thermal resistance is given by: 
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(1.4) 



The overall thermal resistance can be measured experimentally; the wall resistance, 
R^, is known, and the inside resistance can be determined from a chosen correlation. 
Therefore, the outside heat transfer coefficient, h 0 , can be calculated by: 



K - 



— -R A - — 



U 



(1.5) 



* J 



There is one major drawback to this method. As mentioned earlier, when 
laminar flow heat transfer conditions exist within the condenser tube, it leads to poor 
heat transfer performance and therefore represents the dominant thermal resistance 
within the condenser. This inside resistance then becomes the controlling factor in 
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the determination of the outside heat transfer coefficient. Indeed, if the inside 
resistance is a large proportion of the overall resistance, then the accuracy of the 
calculated outside heat transfer coefficient becomes uncertain. Therefore, the 
accuracy of the outside heat transfer coefficient relies heavily on an accurate 
correlation to determine hj. 

An alternative is to reduce the inside controlling resistance by augmenting the 
inside heat-transfer coefficient Augmentation techniques for lowering the inside 
thermal resistance primarily focus on repeatedly bringing fresh fluid from the coolant 
bulk stream to the heated inside surface. This results in breaking up the thermal 
boundary layer and thereby reducing the thermal resistance. The technique utilized 
in-the-NPS-condensation experiments to break up-this thermal boundary layer made 
use of two types of inserts, a twisted tape insert and a wire mesh (HEATEX) insert. 
These inserts and their use will be discussed in more detail in later chapters. 

The results of condensation experiments conducted by Mazzone [Ref. 6] show 
signs of inaccuracies related to the choice and accuracy of the inside heat transfer 
coefficient correlation. For example, in attempting to determine the outside heat 
transfer coefficient for a smooth tube with no insert, Mazzone found that the 
correlation used to determine the inside resistance yielded values that were greater 
than the actual measured overall resistance! Obviously, no sensible values of h 0 could 
be determined. Similar problems existed (but to a lesser extent, due to the greatly 
reduced inside thermal resistance) when using various types of inserts. Mazzone 
concluded that "... laminar flow correlations for twisted tape, HEATEX, and no 
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insert conditions must be developed which have better accuracy than those used ... 
allowing for more accurate determination of outside heat transfer coefficients." 
Therefore, the determination of accurate inside heat transfer coefficients is 
paramount in the final outcome of the refrigeration condensation project 



C. OBJECTIVES 

The main focus of this thesis was therefore to accurately determine inside heat 
transfer coefficients in support of refrigerant condensation experiments being 
conducted at NPS. 

The specific objectives of this study were to: 



1. Design and build a test apparatus to accurately determine inside heat transfer 
coefficients for the same tubes and inserts used in the condensation 
experiments. 

2. Develop correlations for inside heat transfer coefficients under laminar flow 
conditions both with and without inserts. 

3. Using the correlations mentioned above, reprocess existing condensation data 
to determine more accurate values of the outside heat transfer coefficients. 
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II. LITERATURE SURVEY 



A. INTRODUCTION 

Laminar fully developed and developing forced convection flows have been 
analyzed in great detail for various boundary conditions. The study of heat transfer 
in laminar flow through a closed duct was first made by Graetz [Ref. 7] in 1883 and 
later by Nusselt [Ref. 8] in 1910. An excellent and thorough study of the theoretical 
development (with both analytical and experimental solutions) for various laminar 
forced convection internal flow boundary conditions is given by Shah and London 
[Ref. 9]. 

This survey will be concerned with laminar flow heat transfer in fully developed 
and developing flows subjected to both a constant wall temperature and constant 
heat flux boundary condition. It will also summarize some of the augmentation 
techniques used to enhance tubeside laminar flow heat transfer and discuss the 
details involved when using twisted tape and wire mesh inserts to enhance inside heat 
transfer. 
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B. LAMINAR FLOW IN CIRCULAR DUCTS 



1. Boundary Layers and Entrance Lengths 

a. Hydrodynamic Boundary Layer 

For laminar flow in a circular tube, fluid enters the tube with a 
uniform velocity as shown in Figure 2.1; but when the fluid particles make contact 
with the tube wall, viscous effects cause them to assume a zero velocity. These 
viscous effects, associated with shear stresses acting in planes parallel to the tube 
wall, tend to retard the motion of particles in adjoining fluid layers, causing a 
momentum boundary layer to develop. The momentum boundary layer thickness, 5, 
increases until its growth is stopped by symmetry at the tube centerline. The viscous 
effects then extend over the entire cross section of the tube. The flow is then said to 
be hydrodynamically fully developed. At this point, the fluid velocity distribution at 
a given cross section is independent of the axial distance, x (i.e., du/dx=0) [Refs. 9, 
10 , 11 , 12 ]. 

The hydrodynamic entry length for laminar flow (Re D < 2300) is 
defined as that region where the hydrodynamic boundary layer is developing and may 
be obtained from the approximation [Ref. 10]: 






X, 



fdjl 



lam 



0.0 5Re D 



( 2 . 1 ) 
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This expression shows that for Re D =2300, the entrance length is 115 times the tube 
diameter, a considerable length for most experimental facilities. 




Hydrodynamic entrance region region 



Figure 2. 1 Laminar Hydrodynamic Boundary Layer Development in a Circular Tube 
(Courtesy [Ref. 10]) 



b. Thermal Boundary Layer 

When a fluid with a uniform temperature distribution enters a circular 
duct with a wall temperature that is greater than the uniform fluid temperature, 
temperature gradients in the fluid develop and convective heat transfer occurs, as shown 
in Figure 2.2. The region of fluid where these temperature gradients exist is called the 
thermal boundary layer. The thermal boundary layer thickness, 5 t , grows until a 
dimensionless fluid temperature distribution is independent of x. The thermal boundary 
layer is then considered fully developed. This point is usually expressed in 
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Surtact condition 




Figure 2.2 Thermal Boundary Layer Development in a Heated Circular Tube 
(Courtesy of [Ref. 10]) 

terms of a non-dimensional temperature profile given by [Ref. 11]: 



to T,(x)-T m (x) 



( 2 . 2 ) 



where T, is the tube surface temperature, T is the local fluid temperature, and T^ 
is the mean temperature of the fluid over the cross section of the tube. The exact 
shape of the fully developed temperature profile is dependent upon the boundary 
conditions imposed, either a uniform surface temperature or constant heat flux. 

For laminar flow, the thermal entry length may be expressed as the 
following [Ref. 10]: 



D 



0.0 5Re D Pr 



(2.3) 



It should be noted that the rate of development of the velocity and 
temperature profiles in the entry region depend upon the fluid Prandtl number. For 
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Pr = 1, both velocity and temperature profiles develop at the same rate. However, 
for Pr > 1, the hydrodynamic boundary layer develops more rapidly than the thermal 
boundary layer (X fih < X fit ). For very large Prandtl number (Pr > 100), X fih is so 
much smaller than X fit that it is reasonable to assume that a fully developed velocity 
profile exists throughout the thermal entry region. [Refs. 9, 10, 11, 12] 

2. Fully Developed Flow 

Fully developed laminar flow of an incompressible fluid through a circular 
duct exists when both the thermal and velocity profiles are fully developed. Two 
boundary conditions must now be addressed, that of a uniform wall temperature and 
constant wall heat flux. 

For these conditions, the general form of the energy equation in cylindrical 
coordinates for an incompressible fluid having constant properties (with no heat 
generation and neglecting body forces, pressure gradients and viscous dissipation), 
is given by [Ref. 11]: 



It can be shown that for these two cases, the energy equation reduces to 
the following differential equations. 



u — +V iI - — — lr dT ) 



r — 



(2.4) 



clx dr r dr 1 dr 



Constant Wall Heat Flux: 1 d 

r dr 





(2.5) 



a 1 dx 



12 



( 2 . 6 ) 



Uniform Wall Temperature: }_ i r ^£\ - ii 


f T _ T > 


( ^ \ 
*T m 


r dr 1 dr J a 


T r T m 
* m ) 


1 * J 



For the constant wall heat flux boundary condition, Kays [Ref. 11] solves 
the differential equation (equation 2.5) for a fully developed velocity and 
temperature profile and expresses the Nusselt number as: 

Nu = — = 4.364 (2.7) 

k 

Note that for this specified case, the heat transfer coefficient, h, depends 
only on the thermal conductivity, k, and the tube diameter, D. 

A solution to the uniform wall temperature case is more complex, and 
involves evaluating the infinite series solution to the differential equation (equation 
2.6) for the hydrodynamically developed, thermally developing flow situation. In 
evaluating the series solution, the resulting asymptotic Nusselt number is given by 
[Ref. 11]: 

Nu = 3.657 ( 2 - 8 ) 

This result is 16% less than the solution for the constant wall heat flux boundary 
condition. The asymptotic values in equations 2.7 and 2.8 are shown in Figure 2.3. 
Note that the inverse of the Graetz number is equivalent to the dimensionless axial 
distance x‘. 
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3. Entry Length Problem 



The solution to the energy equation for the entry region is far more 
difficult to obtain because velocity and temperature vary both axially and radially. 
There are two entry length problems to be considered, a combined entry length, 
where both temperature and velocity profiles develop simultaneously, and the 
thermal entry length problem, where the thermal conditions develop in the presence 
of a fully developed velocity profile. This latter case would exist if the place at which 
heat transfer begins is preceded by an unheated starting length, a situation that exists 
in most experimental facilities. It should be noted that for Pr << 1 (i.e., liquid 
metals), the thermal boundary layer would develop faster than the hydrodynamic 
boundary layer. However, this case has little practical application and is therefore 
rarely considered in the literature. 

Solutions have been obtained for both entry length conditions mentioned 
above [Refs. 9, 10, 11, 12] and these results are also presented in Figure 2.3. The 
Nusselt number is plotted against a dimensionless axial coordinate x* (the inverse of 
the Graetz number) given by: 



x ' = 



x/D 

Refr 



(2.9) 



Fully developed conditions are reached at x* =* 0.05 where the solutions for the 
thermal entry length converge with those for the combined entry length. For x* > 
0.05, the two types of entry length asymptotically reach the fully developed solutions 
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of 4.36 and 3.66 for the constant wall heat flux and uniform wall temperature 
boundary conditions respectively. 




4.36 

3.66 



Figure 2.3 Local Nusselt Number Obtained from Entry Length Solutions for 
Laminar Flow in a Circular Tube (Courtesy [Ref. 10]) 



a. Thermal Entry Region - Uniform Wall Temperature 

The thermal entry length problem, first studied by Graetz in 1883 and 
later by Nusselt in 1910, is more commonly known as the Graetz Problem. They 
considered an incompressible fluid with constant properties flowing through a 
circular tube having a fully developed laminar velocity profile and a developing 
laminar temperature profile subjected to a uniform wall temperature. Shah and 
London [Ref. 9] present a numerically derived infinite series solution to the Graetz 
Problem. 
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A more useful correlation presented by lncropera and Dewitt [Ref. 
10] for the uniform wall temperature thermal entry length problem was developed 
by Hausen and is of the form: 



0.0 66S(D/L)Re D Pr 

Nu_ = 3.66+ 

U0M[(D/L)Re D Pr] 2/3 



( 2 . 10 ) 



This correlation was used by Mazzone [Ref. 6] in his refrigerant condensation 
experiments to predict the average inside heat transfer coefficient for a smooth tube 
without any type of insert. 

Shah and London [Ref. 9] present a more complete set of correlations 
derived from experimental data for heat exchanger tubes in laminar flow for local 
(Nuj) and mean (Nu m ) Nusselt numbers based on various x* ranges. The correlations 
for local Nusselt number are given by: 



Nu x = \.Qm(x*y V3 -0.1 for x‘ < 0.01 

Nu x = 3.657 +6.874(10 3 x*)-°- 488 e- 57 - 2 *' for x' > 0.01 

and for mean Nusselt numbers by: 



Na m = 1.615(x‘)- l/3 -0.7 
Nu m = 1.615(x*)' 1/3 -0.2 



Nu =3.657 + 



0.0499 



for x ’ ^ 0.005 
for 0.005 <*’ < 0.03 

for x' > 0.03 



( 2 . 12 ) 



The reader is reminded that both the Hausen correlation and the 
Shah and London correlations presented above are for a thermally developing flow 
with a fully developed velocity profile subjected to a uniform wall temperature. 
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b. Thermal Entry Region - Constant Wall Heat Flux 

For the thermal entry length problem under a constant wall heat flux 
boundary condition, Shah and London [Ref. 9] present several numerical solutions 
but recommend the following approximate equations for the local Nusselt number: 

Nu x = 1.302(x*)' u3 -l for x' <, 5X10' 5 

Nu x = 1302Cx*)- w -0J for 5X10' 3 <; x* <; 1.5 xlO* 3 (2.13) 

Nu x = 4.364 +8.68(10 3 x*)' a306 e * 4U ‘ for x' > 1.5X10' 3 

Shah and London also recommend the following approximate equations for the mean 
Nusselt number: 



Nu m = 1.953(x*)' l/3 for x* < 0.03 

0 0722 

Nu m = 4.364 + — — for x' > 0.03 

x' 



(2.14) 



The above two sets of equations are for a thermally developing, 
hydrodynamically developed flow subjected to a constant wall heat flux. 



c. Combined Entry Length 

As discussed earlier, the rate at which the temperature and velocity 
profiles develop is dependent upon the Prandtl number. At medium to high Prandtl 
numbers (Pr > 100), the velocity profile develops much more rapidly than the 
temperature profile and a fully developed velocity profile assumption is valid; for 
such conditions the equations given above should be used. 

For 1 < Pr < 100, however, a fully developed velocity profile is no 
longer a valid assumption and a combined entry length solution must be developed. 
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Again, Shah and London [Ref. 9] present a comprehensive study of various solutions 
to the combined entry length problem, including two solutions for local Nusselt 
numbers developed by Churchill and Ozoe [Ref 13], the first, for a uniform wall 
temperature: 



Nu x + 1.7 

5.357{1 + [(388/tt)x*]- 8/9 P /8 





> 


4/3 


1 + 


7r/(284x‘) 






(1 +(Pr/0.0468) 2/3 ) l/2 {l + [(388/ir)x*]-^P /4 j 





(2.15) 



and the second, for a constant wall heat flux: 



1 + 



Nu x + 1 

5.364(1 +[(220/tt)x*]' 10/9 P /10 
i7/(l 15.2*-) 

{1 +(Pr/0.0207) 2y3 } 1 ' 2 {1 +[(220/ir)x *]- l<w P /3 



5/3 



/ . 



3/10 



(2.16) 



From the preceding paragraphs, it can be seen that the subject of 
laminar flow forced convection in circular ducts has been studied in great detail 
yielding both analytical and experimental results. The reader is referred to Shah and 
London [Ref. 9] and Kays and Crawford [Ref. 11] for a more complete and thorough 
review of this subject. 
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C. AUGMENTATION TECHNIQUES TO ENHANCE TUBESIDE HEAT 
TRANSFER 

1. General Introduction 

As discussed earlier, when laminar flow heat transfer occurs within tubes, 
it usually represents the dominant thermal resistance in a heat exchanger. These 
large resistances result in low inside heat transfer coefficients. In recent years, the 
requirement for more efficient heat exchangers has stimulated considerable interest 
in techniques to augment or enhance tubeside heat transfer. 

The various methods to enhance tubeside heat transfer can be classified 
into two categories, passive and active techniques. Passive techniques involve no 
external energy (except increased pump work) to produce the necessary 
enhancement Passive techniques are divided into the following methods: 



1. Internally finned tubes (longitudinal or spiral fins) designed to increase the 
inside heat transfer surface area. 

2. Surface roughness techniques aimed at agitating flow rather than increasing 
heat transfer surface area. 

3. Swirl flow techniques which involve the use of inlet vortex generators, 
periodically spread propellers or twisted tape inserts designed to cause 
enhancement by creating a rotating and/or secondary flow. 

4. Displaced promoters which alter the fluid flow near the surface by the use of 
inserts such as wire mesh, static mixer elements, rings or disks. 

5. Coolant additives designed to increase the conductivity of the coolant. 

6. Compound techniques that involve more than one of the above mentioned 
methods. 
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Active techniques utilize an external energy source to promote the 
required augmentation. These techniques include: (1) mechanical aids, (2) heated 
surface vibration, (3) fluid pulsation, (4) electrostatic fields, and (5) suction or 
injection. 

All of the above techniques are aimed at creating turbulence in the 
relatively slow moving or stagnant boundary layer fluid by displacing this fluid from 
the heated surface and mixing it with the bulk fluid. Essential in all of these 
techniques is a balance between increased heat transfer rate and increased pumping 
power requirements/costs. Bergles and Joshi [Ref. 14] present an excellent review 
of augmentation techniques and compare heat transfer data taken at low Reynolds 
number flow subjected to both uniforarwaH- temperature and constant wall heat flux. 

Of particular interest to this thesis are two passive techniques, the swirl 
flow technique using twisted tape inserts and the displaced promoter technique using 
wire mesh inserts. These two methods of enhancement of tubeside heat transfer are 
discussed in more detail in the following sections. 

2. Heat Transfer Enhancement with Twisted Tape Elements 

Among the most common swirl flow augmentation techniques is the use 
of twisted tape inserts. The tapes usually consist of a thin stainless steel, brass or 
copper strip which is the width of the inside diameter of the tube. The tapes are 
twisted by clamping one end to an overhead, attaching weights to the other end and 
then twisting to the desired specifications [Ref. 15]. Enhancement occurs not only 
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due to the increased path length of flow, but also due to secondary flow effects and 
fin effects caused by tape contact with the inside tube wall. 



Analytical studies of twisted tape inserts by Date and Singham [Ref. 16] 



suggest that swirl flows improve heat transfer by as much as a factor of 70; however, 
no experimental data have confirmed their unrealistically large predictions. The first 



Bergles [Ref. 15]. Their primary objective was to develop experimentally based 
correlations for predicting heat transfer coefficients for laminar flow of water and 
ethylene glycol in an electrically heated tube (heated length of 1.22 m, inside 
diameter of 10.2 mm) with two sizes of twisted tape inserts. They showed that the 
Nusselt nu mber wa s independent of axial location, .suggesting that fully developed 
conditions existed. Therefore, the mean Nusselt number should take the form: 



experimental work of swirl flow using twisted tapes was conducted by Hong and 




(2.17) 



Re, is the Reynolds number for the swirl flow given by: 




4 m 



(2.18) 



and 5 is the tape thickness. The twist ratio, y, is defined as: 




(2.19) 
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where H is the pitch for a 180° twist of the tape. Hong and Bergles then combined 

Re, and y as a single parameter and came up with the following correlation for 

♦ 

values of Re,/y > 10: 

Nu m = 0.383Pr 0JS (ReJy) 06Z2 (2.20) 

For values of Re/y < 10, the following correlation was obtained: 

Nu m = 5.172(1 +5.484 x 10 ^Pr Q1 {ReJy) l25 ) 0JS (2.21) 

It should be noted that all fluid properties in equations 2.20 and 2.21 
should be evaluated at the bulk fluid temperature. One additional important 
observation was made by Hong and Bergles regarding tube wall temperature. They 
found a significant variation of temperature at one axial location with two different 
tape orientations at that position of 90° and 120°, suggesting that the temperature 
profile was related to the tape orientation. Hong and Bergles concluded that for a 
fully developed flow, inside heat transfer coefficients depend on twist ratio and could 
be improved by a factor of two or three over an empty tube. These conclusions were 
confirmed by later studies of Saha et. al. [Ref. 17] and Mamer and Bergles [Ref. 
18]. 

3. Heat Transfer Enhancement with Wire Mesh (HEATEX) Elements 

The use of wire mesh inserts to promote inside heat transfer falls into the 
displaced promoter technique mentioned earlier. One such wire mesh insert was 
commercially developed by CAL-GAVIN Limited (UK) and called "HEATEX." 
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This HEATEX insert consists of a central wire core on which a series of 



wire loops or petals are attached. The petals are attached to the core at an inclined 
angle that face the oncoming fluid flow. Ideally, each petal makes contact with the 
tube wall in the form of an arc, thereby maximizing disturbance of fluid along the 
tube wall. In addition to the boundary layer disturbance, radial mixing occurs caused 
by repeated redistribution and mixing of fluid through the mesh of loops near the 
tube axis [Ref. 19]. Various wire loop densities and nominal diameters are 
available, depending on the application. 

Initial studies by Gough et al. [Ref. 19] with HEATEX inserts suggest 
enhancements of a factor greater than ten times that of a similar tube with no insert 
Additional experiments were conducted by Oliver and Aldington [Ref. 20] in which 
a correlation was determined: 

Nu = 0 .232Re 0M Pr 0M ( 2 - 22 ) 

It must be remembered, however, that this correlation is for a specific loop density 
and nominal diameter HEATEX element. 

The two types of HEATEX mixing elements supplied by CAL-GAVIN 
(nominal diameters of 13.26 and 10.16 mm) did not match the loop densities of 
previous studies; therefore new correlations needed to be developed. Mazzone [Ref. 
6] developed correlations for these inserts from limited data supplied by CAL- 
GAVIN. For the 13.26 mm and 10.16 mm elements, these correlations were 
respectively: 
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Nu = 0.226 Re 065 Pr 0M 



(13.2 6mm) 



(2.23) 



and 

Nu = 0.0 63 Re 076 Pr° 46 (10.16mm) ( 2 - 24 ) 

Since these correlations were developed from a very limited amount of 
data taken under conditions that did not exactly match that of Mazzone’s 
condensation experiments, they were always subject to inaccuracy. Therefore, the 
need for more accurate correlations under the exact conditions used in the 
condensation experiments are required. 
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III. EXPERIMENTAL APPARATUS 



A. SYSTEM OVERVIEW 

The test apparatus utilized is shown in Figures 3.1 and 3.2. The apparatus is a 
closed loop, low pressure system utilizing an ethylene glycol/water mixture as the 
coolant The coolant is circulated by a centrifugal pump from a supply tank, through 
a reducing manifold and flowmeter into the test section. Coolant is then returned to 
the supply tank via a return header. 

The test section consisted of an inlet and outlet thermocouple well/mixing 
chamber to measure inlet and outlet coolant temperatures and the test tube itself. 
The tube was instrumented with twelve thermocouples (embedded in the tube wall) 
to determine inside wall temperatures. A constant heat flux boundary condition was 
achieved using a series of flexible electrical heating elements wrapped around the 
tube. Monitoring the system as well as data acquisition and data processing were 
provided by a Hewlett-Packard microcomputer and Data Acquisition System. 
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EXPERIMENTAL APPARATUS 
PRIMARY AND SECONDARY COOLANT SYSTEMS 



H 




Figure 3.1 Primary and Secondary Coolant Systems. 
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Figure 3.2 Data Acquisition and Control Systems. 
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B. TUBES AND INSERTS TESTED 



1. Tubes Tested 

Three types of tubes were utilized in the experiments, a copper smooth 
tube, a copper/nickel externally finned tube, and a copper/nickel corrugated or roped 
tube. These tubes correspond to the tubes utilized in the refrigerant condensation 
experiments of Mabrey [Ref. 5] and Mazzone [Ref. 6]. 

The smooth copper tube provided a baseline set of experimental data 
which were then compared with previous experimental and theoretical work for the 
inside heat transfer coefficient conducted by Shah and London [Ref. 9] for no insert, 
Oliver and Aldington [Ref. 20] for the HEATEX insert, and Hong and Bergles [Ref. 
15] for the twisted tape insert 

The Cu/Ni finned tube provided a smaller inside diameter to be used in 
these comparisons. Note that the fins play little part in the experiments. The wall 
thermocouples were buried in the tube wall and were used to measure the inside wall 
temperature, which had a smooth bore. Consequently, for the purposes of these 
experiments, the tube behaved like a smooth copper tube with a different internal 
diameter. 

The final tube tested was a Cu/Ni corrugated (or roped) tube, 
commercially referred to as a "Korodense" tube. The tube was manufactured and 
provided by Wolverine Tube Company of Decatur, Al. All tube dimensions and 
specifications were taken from the Wolverine Technical Bulletin No. 4020 [Ref. 
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21]. Figure 3.3 is a schematic of the Korodense tube. Tube specifications for all 
three tubes are listed in Table E.1 of Appendix E. 




Figure 3.3 Wolverine Korodense Tube (Courtesy [Ref. 21]) 

2. Insert Elements Tested 

Two types of insert mixing elements were utilized in this investigation, a 
twisted tape insert and a wire mesh insert commercially referred to as a "HEATEX" 
tube insert In both cases, two sizes were required due to the two different tube 
inside diameters used. 

The twisted tape inserts were made of a strip of brass 0.559 mm thick and 
either 10.16 mm or 13.26 mm wide. The tapes were manufactured in the NPS 
Mechanical Engineering machine shop by clamping them to the overhead, attaching 
a weight to the opposite end, and twisting them to the desired twist ratio. A twist 
ratio of 3 and 4 corresponding to the two respective tube inside diameters were 
manufactured. Once installed in the tube, the edges of the tape should ideally have 
been in contact with the tube wall. 
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The second type of insert used was a commercially available product 
referred to as HEATEX radial mixing elements. These elements were supplied by 
Cal Gavin of Birmingham, U.K. They are manufactured from stainless steel and 
consist of a central wire core onto which a series of wire loops or petals are attached. 
Each petal was inclined at an angle facing the oncoming flow. The outside diameter 
of the elements were slightly oversized for the tube, ensuring that close contact with 
the tube wall in the form of an arc, rather than a point, was achieved. As with the 
twisted tape elements, different diameter HEATEX elements were needed, 13.26 
mm used in the smooth and Korodense tubes and the other 10.15 mm for the Cu/Ni 
finned tube. Figure 3.4 shows both the twisted tape and HEATEX inserts. 








Figure 3.4 Twisted Tape and HEATEX Inserts 
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C. COOLANT SYSTEM 



1. Primary Coolant System 

An ethylene glycol/water mixture (54/46% by weight) was stored in a 1.5 
m 3 insulated supply tank. The mixture was circulated through the primary system by 
a IV 2 hp, constant speed, centrifugal pump. The coolant was discharged from the 
pump via a 76 mm diameter PVC pipe to a manifold where it was reduced to a 15.9 
mm diameter Tygon flexible tube. A rotameter type flowmeter with a ball valve at 
its entrance was used to control the flow rate through the system. Flowmeter 
calibration results are discussed in Appendix A. 

Flow entered the test section through a 90° bend which housed the 
thermocouple well and exited through a mixing chamber/lhermocouple well 
combination (see Figure 3.1). Prior to entering the exit mixing chamber, the coolant 
passed through a 10 cm length of flexible tubing containing a high density wire mesh. 
This wire mesh, in conjunction with the mixing chamber, served to minimize radial 
temperature variations and provide a uniform outlet coolant temperature. The flow 
then exited the test section through more flexible tubing to a manifold and return 
PVC header, suspended above the apparatus, and back to the coolant supply tank. 

2. Secondary Coolant System 

In order to achieve a wide range of Reynolds number, coolant inlet 
temperatures were varied from -10°C to +20°C. This was accomplished by a 
secondary coolant system comprised of a 28 kW (8 ton) external refrigeration system. 
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A 30 gpm, .75 hp centrifugal pump continuously circulated coolant between the 
supply tank and a chiller barrel where it was cooled to the desired temperature. 
Temperature was maintained by a thermostatically controlled solenoid valve at the 
discharge side of the chiller barrel. This system allowed for a range of coolant 
temperatures between -20° C and ambient to be selected while maintaining the 
coolant inlet temperature approximately constant (±2°C) during each run. 

D. INSTRUMENTATION AND TEST SECTION 
1. Instrumentation 

All temperature measurements were made using copper/constantan type-T 
teflon coated thermocouples and were previously calibrated by Mabrey [Ref. 5]. The 
inlet coolant temperature was measured by a single thermocouple while the outlet 
coolant temperature was measured by two thermocouples located immediately after 
the mixing chamber. 

Each test tube was fitted with twelve wall thermocouples to determine the 
average inside wall temperature. Thermocouples were placed at four different 
longitudinal positions along the axis of the tube. Three thermocouples at each 
location were spread evenly around the tube at 120°C intervals. Figure 3.5 shows a 
typical tube with the exact longitudinal and nodal thermocouple locations. Actual 
longitudinal locations varied slightly for each tube and are given in Table E.l as x/D; 
where x is the length from the start of the heated length to the thermocouple 
location and D; in the inside diameter of the tube. The two extreme longitudinal 
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locations of the thermocouples corresponded to the beginning and end of the 
"effective" condensing length of the tubes being used in the condensation experiments 
of Mabrey [Ref. 5] and Mazzone [Ref. 6]. 



INSTRUMENTED TUBE 



inlet Do 




Di Outlet 



7 




1 1 

POS 1 POS 2 

I5*2cm 55.9cm 



L - 1.52m 
1^= 1.42m 
L/D = 105.6 







5 



POS 3 POS 4 

96.5cm 137.2cm 

Do = 15. 88mm 
Di = 13. 26mm 
Dt/c*14.36mm 



Figure 3.5 Instrumented Tube 



2. Thermocouple Attachment Method 



In order to attach the thermocouples to the wall of the tube, short 
longitudinal grooves were cut at each location. These longitudinal grooves were cut 
to a depth of Vz the tube wall thickness. It was decided that a thermocouple 
embedded within the tube wall vice being attached to the surface would lead to a 
more accurate tube wall temperature and minimize the effect of direct contact 
heating from the heating elements. Thermocouples were then placed in tnese grooves 
and either soldered in place or tack welded and covered with a liquid metal epoxy. 
Different methods for attaching and covering the thermocouples were used for each 
tube based on the tube material and configuration (discussed below). 

In the case of the smooth copper tube, the thermocouples were silver 
soldered in place and the remainder of the groove then covered with silver solder. 
Since the tubes were copper, this was relatively easy to accomplish. Care was taken 
to ensure that the thermocouple remained in contact with the tube wall and that no 
air pockets were created (which would cause an additional large resistance to heat 
flow). Also, great care was taken to prevent any melting of the teflon insulation past 
the thermocouple joint. This condition could have caused an electrical contact 
between the two wires of the thermocouple prior to the tube wall, leading to an 
unrealistically high observed wall temperature. 

A different method was required to attach the thermocouples to the 
copper/nickel finned tube due not only to the external fins but also the 10% nickel 
content. Initially, a small diameter hole was machined into the tube to minimize fin 
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damage (rather than cutting a longitudinal groove). Since the thermocouple was to 
be embedded half way into the tube wall, this hole proved to be relatively deep. 
Also, because of the nickel content of these tubes, the preheat temperature required 
to ensure that the silver solder would adhere to the tube wall was much higher. 
These higher required temperatures and the hole depth caused a significant amount 
of the thermocouple teflon insulation to be melted, allowing for potential contact 
between the thermocouple wires at some point other than the tube wall. This 
technique was therefore abandoned in favor of the longitudinal groove method. 
However, thermocouple insulation was again being destroyed due to the higher 
required preheat temperatures. An alternative solder, soft solder, with lower preheat 
and melting temperatures, was therefore used. The final combination of longitudinal 
grooves and soft solder was successful. 

For the copper/nickel Korodense tube, the groove and soft solder method 
was again tried. Initial tests in the apparatus showed wall temperatures at one 
longitudinal position varying by as much as 15°C. On closer examination of the 
thermocouple connections (using a magnifying glass), it was revealed that insulation 
damage was evident on over half of the thermocouples. Several more attempts at the 
groove and soft solder method were tried, all leading to the same results. 

A third method utilizing a tack welding machine (acquired and adapted 
for our use) was tried. Thermocouples were tack welded into the bottom of the 
groove resulting in no damage to the insulation and a liquid steel epoxy was used to 
fill the remainder of the groove, providing additional strength and protection to the 
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connection. This final method proved to be the most successful and is recommended 
for use in any follow-up experiments. 

3. Test Section 

In order to achieve a constant heat flux, each tube was tightly wrapped 
with six silicone rubber flexible electrical heating elements rated at 240 V, 200 W 
each. Each strip was 2.5 cm wide and 50.6 cm long. The heating elements were 
connected in parallel to a 240 Vac power supply via a plug panel and power 
switchboard. Power to the six heaters was controlled by a STACO, 240 Vac, 23.5 kVa 
rheostat (see Figure 3.2). 

To minimize heat loss from the heater to ambient, the tube was wrapped 
in a 20 mm thick neoprene type insulation. This type of insulation proved to be 
extremely effective giving a negligible heat loss. It also provided insulation against 
heat influx from ambient at the colder coolant inlet temperatures. Actual heat loss 
and heat influx experiments and their results are discussed in more detail in 
Appendices B and C. All other tubing and PVC piping was wrapped in similar type 
insulation. 

The whole test section was then placed in a restraining stand which 
minimized outside vibration effects and ensured that each tube remained horizontal 
and in the same relative position (to the remainder of the apparatus) for each 
experiment. 
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E. DATA ACQUISITION AND CONTROL SYSTEM 

A Hewlett-Packard 9125 microcomputer (MC) was used to control a Hewlett- 
Packard 3497A Automatic Data Acquisition Unit (DAU) (see Figure 3.2). The DAU 
read the output from the twelve tube wall thermocouples and three coolant 
thermocouples (one inlet and two outlet). Additional channels were assigned to sense 
heater input voltage and current at the plug panel. The voltage was converted to 
heater power using the heater resistance and the calibration discussed in Appendix 
D. 

Thermocouple readings were made in microvolts which were converted to 
temperatures (along with the heat flux calculations) using the data reduction 
program. The data reduction program, DRPSING, will be discussed in more detail 
in a later chapter. All temperatures, power, and heat flux measurements were 
monitored using the MC and DAU. DAU channel assignments are listed in Table 
3.1. 
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TABLE 3.1 DATA ACQUISITION UNIT CHANNEL ASSIGNMENTS 



CHANNEL NO. 


MEASURED QUANTITY 


1 


Thermocouple Pos 1 


2 


Thermocouple Pos 2 


3 


Thermocouple Pos 3 


4 


Thermocouple Pos 4 


5 


Thermocouple Pos 5 


6 


Thermocouple Pos 6 


7 


Thermocouple Pos 7 


' 8 


Thermocouple Pos 8 


9 


Thermocouple Pos 9 


10 


Thermocouple Pos 10 


11 


Thermocouple Pos 11 


12 


Thermocouple Pos 12 


13 


Tube Inlet 


14, 15 


Tube Outlet 


25 


Input Current 


29 


Input Voltage 
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IV. EXPERIMENTAL PROCEDURES 



A. CALIBRATION EXPERIMENTS 

Several calibration experiments were required prior to performing the main 
experiments. Figure 4.1 depicts a typical tube being tested. The primary heat source 
is from the heating elements (q up< ). However, to complete the heat balance on the 
system, several other heat sources and heat sinks needed to be accounted for. The 
first experiment conducted was to determine the amount of heat being lost from the 
heating elements to the atmosphere vice going to the coolant Heat loss experiments 
are discussed in Appendix B. 

The second experiment was to account for two additional potential heat 
sources: heat leakage into the coolant from ambient (q^) and frictional temperature 
rise (q fncl ). Approximately 8 cm at either end of the tube, along with the flexible 
tubing section prior to the thermocouple wells, were not covered by the heating 
elements. Since coolant inlet temperatures were primarily below ambient 
temperature, a potential for heat leakage to the coolant from ambient existed. Also, 
because of these cold inlet temperatures, frictional heating could exist, especially 
when using the various inserts. Heat leakage into the system and frictional effects 
were considered together as one heat influx. These experiments were conducted with 
the heater tapes turned off; any temperature rise was then monitored using the inlet 
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and outlet thermocouples. Appendix C discusses the experiments conducted to 
determine this part of the heat balance equation. As mentioned earlier, the test 
section and all other tubing were insulated to minimize these effects. 



Figure 4.1 Heat Balance. 

B. PROCEDURES FOR DETERMINATION OF INSIDE HEAT TRANSFER 
RATE 

1. General Procedures 

Each tube utilized in this study was tested with no insert and the 
appropriate twisted tape and HEATEX inserts corresponding to the tube inside 
diameter. For each insert condition, four nominal coolant inlet temperatures, -10°C, 
0°C, + 10°C and +20°C were chosen to increase the range of experimental Reynolds 
and Prandlt numbers that could be attained. To check that heat flux variations did 
not affect the results, data were taken at heat fluxes of 10 kW/m 2 and 15 kW/m 2 for 
each inlet temperature. For each insert, coolant inlet temperature and heat flux 
setting, coolant velocity was varied between 0.2 m/s and 1.4 m/s. Since the inside 
diameter of the Cu/Ni finned tube was smaller than the other two tubes, higher 
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coolant velocities of up to 2.4 m/s could be attained. A matrix of the individual data 
runs and data file names is shown in Table 4.1. 

With the tube installed and all calibration experiments completed, the 
experiments to determine the inside heat transfer coefficient were conducted. Since 
flowmeter settings for specific coolant velocities were based upon the inlet coolant 
temperature, the system was initially operated at a 100% flowmeter setting with no 
power to the heating elements for 15 minutes. This allowtd the inlet thermocouple 
well to stabilize. Once the desired inlet coolant temperature was reached and 
maintained, power to the heating elements was set to the predetermined heat flux 
setting for that data run (10 or 15 kW/m 2 ). The flowmeter was then adjusted to 
achieve the desired coolant velocity based upon flowmeter calibration results 
(discussed in Appendix A). Flow was maintained for approximately 10 minutes to 
ensure steady state equilibrium conditions prior to taking data. 



TABLE 4.1 DATA FILE NAMES 



TEMP(C)/ 


Cu Siooth Tube 


Cu/Nl Finned 


Tube 


Cu/Ni Korodense Tube 


q* (kw/» J ) 


NO 

Insert 


Twisted 

Tape 


HEATEX 


NO 

Insert 


Twisted 

Taoe 


HEATEX 


No 

Insert 


Twisted 

Tape 


HEATEX 


20/10 


SUN 1 04 


- 


SWHX04 


CFNI04 
( + 15) 


CFT01 

(+15) 


CFHX01 

(♦151 


KDNI01 


KDTT01 


KDHX01 


20/15. 


SMNI04 


- 


SJ4HX04 


CFNI04 

(-15) 


CFTT01 

(♦15) 


CFHX01 

(♦15) 


K0NI01 


KDTT01 


KDHX01 


10/10 


SMNI03 


SMTT01 


SMHX01 


CFNI01 


CFTT02 


CFHX02 


KDNI02 


KDTT02 


KDHX02 


10/15 


SMNI03 


SMTT01 


SMHX01 


. 


CFTT02 


CFHX02 


KDNI02 


KDTT02 


KDHX02 


0/10 


SMNI02 


SMTT02 


3MHX02 


CFNI03 


CFTT03 


CFHX03 


KDNI03 


KDTT03 


KDHX03 


0/15 


SMNI02 


SMTT02 


SMHX02 


CFNI03 


CFTT03 


CFHX03 


K0NI03 


K0TT03 


KDHX03 


-10/10 


SMNI01 


SMTT03 


SMHX03 


CFNI02 


CFTT04 


CFHX04 


KDNI04 


K0TT04 


K0HX04 


-10/15 


SMNI01 


SMTT03 


SMHX03 


CFNI02 


CFTT04 


CFHX04 


KDNI04 


KDTT04 


KDHX04 



41 



Velocities were increased from 0.2 m/s to 1.4 m/s in steps of 0.2 m/s and 
then decreased from 1.3 m/s to 0.3 m/s in similar steps. This method of increasing 
on even values and decreasing on odd values was done not only to determine if any 
hysteresis effects existed (i.e., a difference in the increasing and decreasing data), but 
also to ultimately increase the number of data points acquired. Also, since data were 
immediately plotted (as Nu m versus Re), the method allowed for an easy 
determination as to whether a data point had reached a steady state condition. 

Inlet and outlet coolant temperatures, tube wall temperatures, heat flux 
and heating element power settings were all continuously monitored using the 
computer program DRPSING in conjunction with the Data Acquisition System 
discussed earlier (DRPSING is listed in Appendix G). Individual data runs for each 
tube and insert condition were conducted in a similar manner. 

2. Specific Experiments 
a. Copper Smooth Tube 

The first tube tested was the copper smooth tube. This tube was 
chosen to provide a baseline set of data needed for validation of system performance 
and for comparison with data acquired using other tubes. It should also be noted that 
two complete independent sets of data covering all insert conditions and inlet 
temperatures were taken over a period of two weeks. Data from both sets were 
compared and found to be in very good agreement (±2%); this validated the 
repeatability of data obtained from the system under the same conditions. 
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The initial intention for all tubes was to start with a given insert at 
the highest coolant inlet temperature, take data at each of the heat flux settings and 
then proceed to the next coolant temperature. Once one insert was completed, the 
system could then be shut down and the insert changed. This procedure would then 
be repeated for all three insert conditions. However, the coolant supply tank and 
coolant system were used concurrently for other experiments which required the 
coolant temperature to be -20°C. Since it took approximately 36-48 hours for the 
coolant temperature to increase from -10°C to +20°C (the sump was not fitted with 
a heater), only the copper smooth tube was done in this way. The method for the 
other two tubes is discussed below. 

b. Cu/Ni Externally Finned Tube 

The second tube tested was the Cu/Ni finned tube. For a given inlet 
temperature and tube insert, the same procedure as with the copper smooth tube was 
used. However, rather than decreasing coolant temperature with the same insert, the 
coolant temperature was maintained and the insert changed until all three of the 
insert conditions were completed for a given inlet temperature. Once all the data for 
a given inlet temperature had been taken, the supply tank temperature was 
decreased to the next inlet temperature condition. It took only about one hour for. 
the secondary coolant system to decrease the temperature by 10°C. Since it was a 
relatively short task to change inserts, the 36 to 48 hour waiting period required to 
raise coolant temperatures back to +20°C between insert changes was eliminated. 
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It should be noted that for the Cu/Ni finned tube, the highest nominal 
inlet temperature attained was only +15°C vice +20°C. Colder than normal ambient 
temperatures existed preventing supply tank temperatures from reaching +20°C. 
Also, due to the smaller inside diameter, coolant velocities could be varied between 
0.3 m/s and 2.4 m/s. 

c. Cu/Ni Korodense Tube 

The same procedure used with the finned tube was also used with the 
Korodense tube. In addition, another test was conducted with this tube. It was noted 
that larger than normal wall temperature variations existed at one longitudinal 
location when using the two inserts. It was suspected that the point of contact 
between the insert and tube wall w as direc tly below a thermocouple position causing 
a higher wall temperature at that point To verify this, the system was shut down by 
securing the coolant pump and the heating elements so that the flowmeter and 
power settings remained exactly the same. The outlet end of the tube was opened 
and the insert moved downstream by approximately 2 cm. The system was then 
closed and the coolant flow and power to the heating elements restored. Wall 
temperature variations still existed but were significantly less. This procedure was 
repeated several times with similar results. 



44 



V. DATA REDUCTION 



A. AVERAGE INSIDE HEAT TRANSFER COEFFICIENT 

The main objective of the data reduction scheme was to determine the average 
inside heat transfer coefficient (hj). This in turn was used to determine a mean 
Nusselt number (for the whole heated length of the tube) to be used as the basis for 
establishing the desired correlations. 

To accomplish this task, the Data Acquisition Unit discussed in Chapter III 
made use of the computer program DRPSING. A listing of the program is given in 
Appendix E. The program allowed the user to monitor information from the 
instrumentation attached to the apparatus, acquire and then process the data to 
determine the desired output In doing this, the user would select the tube and insert 
condition as well as the nominal coolant inlet temperature for that particular data 
run. The desired coolant velocity was then selected. Based on the actual measured 
coolant inlet temperature, coolant mixture properties and flowmeter calibration 
results (see Appendix A), the appropriate flowmeter setting (as a %) was calculated 
and set on the flowmeter. 

The program then monitored coolant inlet and outlet temperatures, the twelve 
wall thermocouple temperatures and the actual measured heat flux and power 
settings. Once all readings reached a steady state condition (approximately 10 min.). 
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the user could terminate the display session and acquire data for that coolant velocity 
setting. 

The DAU measures emf from the thermocouples and then converts these 
values to temperatures using the temperature calibration equation developed by 
Mabrey [Ref. 5]. Heater power was determined by measuring the voltage at the plug 
panel and calculating the power based on heater resistance and the power calibration 
contained in Appendix D. 

The twelve wall thermocouples (placed axially and radially) were used to 
minimize the uncertainty in determining the average inside wall temperature. These 
twelve wall temperatures were averaged and corrected (to account for the depth of 
burial using a simple conduction calculation) to determine the average inside wall 
temperature, T^. 

The physical and thermodynamic properties of the ethylene glycol/water 
mixture (density, dynamic viscosity, thermal conductivity, heat capacity) were 
determined using equations and figures from Cragoe [Ref. 22] and Gallant [Ref. 
23] and correlated by Mabrey [Ref. 5]. From these quantities, the coolant mass 
flow rate was calculated by: 

m = V c pA i (5.1) 

where A ; is the inside cross sectional area of the tube and V c is the set coolant 
velocity. The coolant Reynolds number, Re 0 was calculated by: 
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(5.2) 



Re D = 



v Ap< 



where D; is the inside diameter of the tube. The dimensionless distance (x*) in the 
flow direction is calculated as: 



x ' = 



x/D. 

~Re^Pr 



(5.3) 



where x is the axial distance from the beginning of the heated length. 

To determine the heat transfer rate to the coolant, several factors were 
considered. The basic equation to determine the heat transfer rate is given by: 



£{ - mcJ^T (5.4) 

where AT is the measured coolant temperature rise along the tube. The temperature 
rise required is that for the actual heated length of the tube that is subjected to the 
constant heat flux condition, AT HL . To determine AT HL from AT, corrections have to 
be made to account for heat loss, heat leakage and heat transferred to the unheated 
ends. Figure 5.1 depicts a typical tube and the correction terms considered in the 
heat balance. 

The overall AT is first corrected for heat loss to the ambient, (q losJ ) and heat 
leakage into the system, (q leak ). The heat leakage term consisted of any temperature 
rise due to the influx of heat from ambient as well as that caused by frictional 
effects. Appendix D discusses the experiments conducted to determine both this heat 
loss and heat leakage correction. 
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TEMPERATURE RISE CORRECTIONS 




^tuo* 

* T c 



T 

t/c.out 



Figure 5.1 Temperature Rise Corrections. 

To take account of the heat transferred to the coolant from the unheated ends, 
the heat transfer rate to the coolant over the entire length of the tube (heated and 
unheated) is first calculated. The average heat transfer coefficient can then be 
determined by: 




Qhl 

A~LMTD 



(5.5) 



where Aj is the inside surface area based on the heated length (L H ) and LMTD is 
the log mean temperature difference defined as: 

At this point in the calculation, however, the heat transfer rate to the coolant 
from the heated length alone (q HL ) is not known. Since the heating elements did not 
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LMTD = 



T -T 

<mi * n 


In 


1 

i 

i 




i 

r 1 

1 

r 

i 





(5.6) 



cover the entire length of the tube, account must be taken of heat added to the 
coolant from the unheated lengths at the inlet and outlet of the tube. This was 
accomplished by considering these ends as annular fins. An adiabatic fin tip and 
outer wall boundary condition was applied and the heat transfer from the fins, q f (as 
given by Incopera and Dewitt [Ref. 10]), determined from: 



q f = ^7iPkA e Q b \anhmL f 



where L* = fin length 
P = fin perimeter = irDi 
A e = fin cross sectional area 

4 



0 b = T W i-Tin (for inlet fin) 
= T^-Tou, (for outlet fin) 



m = 
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The fin calculation requires a known average inside heat transfer coefficient (h ; ). 
Therefore, an initial estimate of h ; is calculated by: 




4 tub* 



(5.8) 



A.LMTD 



where q, ub< is known. The amount of heat that passes through the unheated ends is 
then calculated using the fin correction above. The heat transfer rate to the coolant 
from the heated length (c^l) is then calculated by: 



A new average inside heat transfer coefficient (h;) can then be calculated using 
equation 5.5 and an iteration process carried out until convergence to a final average 
inside heat transfer coefficient is reached. 

An additional correction for the Cu/Ni finned tube is needed to account for the 
change in inside diameter of the smooth ends (needed to fit through the tube sheet 
in the condenser). The heat transfer coefficient based on the smaller inside diameter 
corresponding to the finned portion of the tube is what ultimately is desired. As 
shown in Figure 5. 2, the heated length extended beyond the smaller diameter section. 
To correct for this, the annular fin correction discussed earlier was applied to 
determine the temperatures at the inlet and outlet of the heated length. The coolant 
temperature rise along the heated length was assumed to be linear, and therefore the 
temperature rise along the smooth larger inside diameter lengths (L t and L,) could 
be determined by the simple ratio: 



QhL ~ 4tub*~4f 



(5.9) 
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(5.10) 



LAT hl , _ LAT hl 

A 7\ = 1 — ami AT 2 = 2 

where Lhl is the heated length of the tube and AThl, the temperature rise along the 
heated length. These corrections were then added or subtracted from the corrected 
finned inlet and outlet temperatures. 



Cu/Ni EXTERNALLY FINNED TUBE 
TEMPERATURE RISE CORRECTIONS 




inlet *“ AT,— 



AT — outlet 
2 



* L 

i 



Lhl 



L 2 



Figure 5.2 Cu/Ni Finned Tube Temperature Rise Corrections. 

With the average inside heat transfer coefficient known, the mean Nusselt 
number can be determined by: 



Nu. 

fff 




( 5 . 11 ) 



where k is the thermal conductivity of the ethylene glycol/water mixture. 
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The experimentally determined mean Nusselt number was then compared to 
various theoretical and experimental correlations, based on the insert condition. 
These correlations were discussed earlier but are briefly mentioned here. 

For the no insert condition, the data were compared to two theoretical 
correlations, the Hausen correlation [Ref. 10] given by: 



Nu 

n% 



0.0668 (D/L)Re D Pr 
l+0.04[(D/L)Re D Pr] 2 ' 3 



(5.12) 



and the Shah and London correlation [Ref. 9] given by: 

Nu m = 1.93(jc*)' u3 (5-13) 

The reader is reminded that the Hausen correlation is for a uniform wall 
temperature boundary condition and was used in Mazzone’s condensation 
experiments [Ref. 6]. The Shah and London correlation is for a constant wall heat 
flux condition. It should also be noted that equation 5.13 is used since x* for the test 
section was always less than 0.03. 

For the twisted tape insert the experimental mean Nusselt number was 
compared to the correlation developed by Hong and Bergles [Ref. 15]: 



Nu m = 0.383 Pr 035 (ReJy)°- 6n 



(5.14) 



Equation 5.14 is used because values for Re^y were always greater than 10. 



52 



For the HEATEX insert, the data were compared to the two correlations 
developed by Mazzone [Ref. 6]. For the 13.26 mm inside diameter tube, the equation 
is given by: 

Nu m = 0.226 Re 065 Pr 0M (5.16) 

and for the 10.16 mm inside diameter tube by: 

Nu = 0.063 Re 0J6 p r 0M (5.16) 

fit 

B. ADDITIONAL TUBE LENGTHS 

In addition to the above calculations for the entire tube length, mean Nusselt 
numbers for three additional tube "lengths" could be determined. This was 
accomplished by sectioning the heated length of the tube into four cells as shown in 
Figure 5.3. Each cell contained one longitudinal thermocouple position. Three 
additional tube lengths were then determined by adding cell lengths together. The 
first length, L t , consisted of cell 1 alone, Lj consisted of cells 1 and 2 and L, 
consisted of cells 1, 2 and 3. The fourth length was actually the entire heated length 
as previously discussed. 

The coolant temperature rise along the heated length of the tube was assumed 
to be linear. The correction to the wall heat flux to account for variation in electrical 
resistance with temperatures was assumed to be negligible. Therefore, the 
temperature rise along each cell could be determined by: 
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SECTIONED TUBE 





^HL 

T/C POS 1 T/C POS 2 T/C POS 3 T/C POS 4 






INLET 


CELL 1 


CELL 2 


CELL 3 

i 

i 
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i 

CELL, 4 


OUTLET 










Figure 5.3 Additional Tube Lengths. 

- (5.17) 

Average inside wall temperatures were based on the thermocouples within that tube 
"length." Therefore, the temperature rise, heat flux and average inside wall 
temperatures were all known for each tube length. Applying the same calculation 
method as before, the mean Nusselt number for these three additional tube "lengths" 
could be calculated. 

This adaption to the program was not done until after the smooth tube 
experiments had been completed; however, this extra data were taken for the 
remaining two tubes. 
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C. LOCAL NUSSELT NUMBERS 



The program was also adapted to determine local inside heat transfer 

♦ 

coefficients at the four axial thermocouple positions. Corrections discussed earlier 
to determine the temperatures at the inlet and outlet of the heated length were first 
accounted for. Exact thermocouple longitudinal positions were entered into the 
program as L/D ratios so that a local x* could also be calculated. Again, the 
temperature rise along the heated length was assumed to be linear; therefore, the 
bulk coolant temperature at the particular thermocouple position could be 
determined by the simple ratio method discussed earlier (equation 5.10). 

The local heat transfer coefficient, h*, was then determined by: 



h . = 



t ,-r 



(5.13) 



where T^ is the average inside wall temperature based on the three thermocouples 
at that position and T c is the bulk coolant temperature at that position. The local 
Nusselt number was then given by: 




(5.19) 



The local Nusselt number for all four thermocouple positions was then determined. 

For the smooth tube experiments, local Nusselt numbers were calculated only 
for position 1. The program was later adapted to include local calculations at 
positions 2, 3, and 4 for the Cu/Ni finned tube and Korodense tube. 
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VI. RESULTS AND DISCUSSION 



A. GENERAL 

The experimental heat transfer results are initially presented in terms of the 
mean Nusselt and Reynolds numbers for all insert conditions. From these graphs, it 
can be seen how the Nusselt number is inlet temperature dependent, indicating the 
need for a Prandtl number (raised to an exponent) in the final correlation. However, 
although not specifically shown, the data shows no dependence on the heat flux 
setting. In addition, the Nu vs. Re graphs for the no insert condition show when a 

* 

transition to turbulent flow occurs. 

The heat transfer results are then compared to existing correlations based upon 
the insert condition. For the no insert condition, the Nusselt number is plotted 
against x* (dimensionless x coordinate). Note that for a given tube, L/D is constant 
and x* varies due to Re and Pr changes. The data are then compared to the Hausen 
[Ref. 10] correlation (equation 2.10) and the Shah and London [Ref. 9] correlation 
(equation 2.14). Recall that the Hausen correlation is for a uniform wall temperature 
boundary condition and was used by Mazzone [Ref. 6] in the refrigerant 
condensation experiments. Since the present data were taken under a constant heat 
flux boundary condition, comparison with the Shah and London correlation was also 
required. The twisted tape insert data are compared to the correlation developed by 
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Hong and Bergles [Ref. 15] (equation 2.20) which was also used by Mazzone. The 
HEATEX insert data are compared to the two correlations developed from the 
limi ted data received from CAL-GAVIN by Mazzone (equations 2.27 and 2.28). 

Correlations are developed for each tube and insert condition based upon a 
least squares linear regression. An r-squared value, a measure of the accuracy of the 
least squares curve fitting to the data, is used to determine the uncertainty of the 
correlation. The correlations and the corresponding r-squared values are summarized 
in Table 6.1. Note that the 0.2 m/s and 0.3 m/s data were eliminated prior to 
developing all correlations due to the high degree of uncertainty in these data points. 
TABLE 6.1 SUMMARY OF CORRELATIONS 

SUMMARY OF CORRELATIONS 





NO INSERT 


TWISTED 

TAPE 


HEATEX 




Nu « 


Nu * 


Nu ■ 


SMOOTH 


4.l37Xstar~-0.234 


.5l2Pr~35(Res/yr.508 


.104Pr'\40Re~.64l 


TUBE 


r-sq ■ 71.9% 


r-sq ■ 99.5% 


r-sq * 98.9% 


EXT 


Nu • 


Nu ■ 


Nu » 


FINNED 


1.84Xstar~-0.34i 


.148Pr“35(Res/y)‘087 


.l20Pr'.46Re~057 


TUBE 


r-sq * 93.0% 


r-sq * 90.7% 


r-sq * 95.9 




Nu • 


Nu • 


Nu ■ 


KORODENSE 

TUBE 


2.O47Xstar~-O.320 


.509Pr'.35(Res/y)*498 


.282Pr'.46Rg'.499 


r-sq * 39.2% 


r-sq * 98.9% 


r-sq * 94.3% 
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B. SMOOTH TUBE 



1. Smooth Tube • No Insert 

The first tube tested was the copper smooth tube with no insert. Heat 
transfer data is presented in Figure 6.1 in terms of the Nusselt number as a function 
of Reynolds number for each nominal inlet temperature. The first thing to notice is 
the sudden increase in Nu at a Re number of about 2300. This is indicative of a 
change from laminar to transition flow where the heat transfer is greatly increased 
due to improved mixing. This contradicts Mazzone’s [Ref. 6] assumption of a 
transition to turbulence occurring at Re > 4000. The data in the transition region 
was primarily for the +20°C and +10°C nominal inlet temperatures. Figure 6.1 also 
shows distinct groupings of data at given inlet temperatures, indicating a Prandtl 
number effect. These groupings are shown more clearly in Figure 6.2 which represent 
data under laminar flow conditions only (Re ^ 2300) and Figure 6.3 for data in the 
transition region only (Re > 2300). 

To develop a useful heat transfer correlation, the mean Nusselt number 
is plotted as a function of x* in Figure 6.4. Recall that x* can be determined by: 



x * = 



L/D 

Re D Pr 



(2.9) 



For a given tube, LTD is a fixed quantity. Therefore, for a given Prandtl number 
(essentially fixed inlet temperature), x* is dependent only upon the Reynolds number. 
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Figure 6.4 also shows a comparison of the experimental data with the 
Hausen correlation (equation 2.10) for a uniform wall temperature (used by 
Mazzone) and the Shah and London correlation (equation 2.14) for a constant heat 
flux. The Hausen correlation appears to significantly underestimate the heat transfer 
coefficient found in the present experimental data, i.e., the present experimental data 
indicate a lower inside thermal resistance. This is consistent with Mazzone’s [Ref. 6] 
results in that the Hausen correlation was predicting larger inside resistances than 
the overall measured resistance. However, it must always be remembered that the 
Hausen correlation is for a uniform wall temperature boundary condition and this 
would indeed lead to lower heat-transfer coefficients. 

The data also fall slightly above (but parallel to) that predicted by the 
Shah and London correlation. Due to the constant heat flux boundary condition used 
by Shah and London, one may expect better agreement (as seen). However, the fact 
that the data fall slightly above can be attributed to secondary flow induced by the 
apparatus inlet conditions (90° bend due to the thermocouple well). Also, as before, 
a departure from laminar flow can be seen in Figure 6.4. Since the transition data 
forms two distinct lines (+20°C and +10°C), it indicates that x’ is not a good 
correlating parameter for transitional flow. 

After removal of the transitional flow data, a correlation was developed 
based upon the form of the Shah and London correlation. The resulting present 
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correlation for the mean Nusselt number for laminar flow with no insert is given by: 



Nu m = 4.137(jc*)"° -234 (6.1) 

Equation 6.1 and the Hausen correlation, equation 2.10, are both shown 
in Figure 6.5 with the present laminar flow data. The general form of Figure 6.5 is 
explained by the fact that at low x* (high Re) the heat transfer is improved due to 
slightly better mixing. As x* increases (coolant velocity decreases), the heat transfer 
decreases. 

Figure 6.5 also shows a high degree of scatter in the data. Initially, it was 
thought that this scatter was due to the flow not being hydrodynamically fully 
developed. However, it is seen to be more pronounced at large x*. Recall that for a 
given nominal inlet temperature, large x* represents a low Reynolds number. For a 
low Reynolds number, the hydrodynamic entry length (given by X/D = 0.05Re) is 
the shortest. Also, since we are dealing with relatively large Prandtl numbers (30- 
140), the assumption that the velocity profile is fully developed throughout the 
thermal entry region should be fairly valid. Therefore, it is felt that this scatter is not 
due to the flow being hydrodynamically undeveloped. 

This scatter is believed to be due to the effects of the different inlet 
temperatures, indicating that x' is not an ideal correlating parameter. That is to say, 
the Prandtl number exponent needs to be adjusted to account for the different inlet 
temperatures. It should also be noted that due to the high degree of detail in Figure 
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6.5, the scatter appears more pronounced. However, in Figure 6.4 with larger scale 
increments, the scatter seems minimal. 

2. Smooth Tube • Twisted Tape 

Results for the smooth tube with a twisted tape insert are shown in 
Figures 6.6 and 6.7. A definite Praudtl number dependence can be seen in Figure 6.6 
where the mean Nusselt number is presented as a function of Reynolds number. 
Following the guidance of Hong and Bergles [Ref. 15] the dimensionless parameters 
Re, and y were combined as one parameter Rejy and the mean Nusselt number and 
Prandtl number were combined as Nu/Pr OJS . The final correlation shown in Figure 
6.7 is: 



Nu m = 0 JSllPr^iReJyf** ( 6 - 2 ) 

This correlation is in close agreement with the Hong and Bergles correlation 
(equation 2.20). The slight differences can be attributed to differences in the 
experimental apparatus (including the insert itself) and the possible existence of 
secondary flows. Figure 6.7 shows little inlet temperature dependence, indicating that 
the Prandtl number exponent of 0.35 is probably fairly accurate. Furthermore, there 
is little experimental data scatter, giving increased confidence in this correlation. 

3. Smooth Tube - HEATEX Insert 

HEATEX insert results for a smooth tube are shown in Figures 6.8 and 
6.9. Again, on a Nu vs. Re basis, it can be seen from Figure 6.8 that there is a 
definite dependence on the inlet temperature. Using the same form as the 
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correlation of Oliver and Aldington [Ref. 20], the resulting correlation for the mean 
Nusselt number is: 

Nu = 0. 164 Pr 046 Re 0Ml (6.3) 

n% 

This is compared in Figure 6.9 to the correlation developed and used by Mazzone 
(equation 2.23). It is clearly evident that Mazzone’s correlation underpredicts the 
inside thermal resistance and therefore gives a larger inside heat transfer coefficient 
Again, Figure 6.9 shows little dependence on inlet temperature indicating that the 
Prandtl number exponent is fairly accurate. Comparisons of the enhancement in heat 
transfer gained by the use of HEATEX and twisted tape insert, over the no insert 
condition are given section F of this chapter. 
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SMOOTH TUBE - NO INSERT 



Nu vs Re 




a +20C ♦ +10C A OC . -IOC 

SMNIRENU 



Figure 6.1 Smooth tube - No insert (Nu m vs. Re) results. 
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Figure 6.2 Smooth tube - No insert (Nu m vs. Re) laminar flow results. 



2500 



64 



SMOOTH TUBE - NO INSERT 

TRANSITIONAL FLOW 

70 



60 - 



b 

o 



50 - 



o> 

VI 



2 40 - 



c 

<8 

<D 



30 - 



20 L 



s 



10 

2000 



i 

2500 



Nu vs Re 

SMNINURETRAN 



3000 3500 4Q00 

Reynolds 

a +20C ♦ +10C 



4500 



Figure 6.3 Smooth tube - No insert (Nu m vs. Re) transitional flow results. 
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Figure 6.4 Smooth tube - No insert (Nu m vs. x*) results. 
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Figure 6.5 Smooth tube - No insert laminar flow correlation results. 
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Figure 6.6 Smooth tube - Twisted tape (Nu m vs. Re) results. 
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Figure 6.7 Smooth tube - Twisted tape correlation results. 
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Figure 6.8 Smooth tube - HEATEX (Nu^ vs. Re) results. 
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Figure 6.9 Smooth tube - HEATEX correlation results. 
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C. CU/NI EXTERNALLY FINNED TUBE 



1. Cu/Ni Finned Tube - No Insert 

The second tube tested was the Cu/Ni externally finned tube. Recall that 
this tube has a smaller inside smooth diameter than the smooth tube (10.15 mm vs. 
13.26 mm). The no insert results for the mean Nusselt number as a function of 
Reynolds number are shown in Figure 6.10 and in more detail in Figures 6.11 and 
6.12 for laminar and transitional flows, respectively. It can be seen that a slightly 
higher Reynolds number (Re = 2700) for the departure from laminar flow exists for 
this tube. Again, the dependence on inlet temperature is clearly evident 

The data are compared to the Hausen correlation (equation 2.10) and the 
Shah and London correlation (equation 2.14) in Figure 6.13. As expected, the 
Hausen correlation again underpredicts the inside heat transfer coefficient. However, 
the laminar flow data is in closer agreement with the Shah and London correlation 
than they were for the smooth tube. This can be attributed to the different diameter 
affecting the thermal boundary layer development. The resulting correlation for the 
laminar flow mean Nusselt number is shown in Figure 6.14 and is given by: 

Nu n = 1.84 (x*)' 0 - 541 (6.4) 

The same conclusions about data scatter can be drawn as for the smooth tube 
results, although in general the scatter seems to be less. 
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2. Cu/Ni Finned Tube - Twisted Tape Insert 

Figures 6.15 and 6.16 show the results for the Cu/Ni finned tube using a 
twisted tape to enhance tubeside heat transfer. The inlet temperature dependence 
in Figure 6.15 seems to be reduced, but still evident The form of the correlation was 
again that of Hong and Bergles’ correlation. It was developed in the same manner 
as for the smooth tube and is given by: 

Nu m = 0. 148Pr 0J3 {Re Jy) 0 - 687 (6.5) 

The data and resulting correlation are compared to the Hong and Bergles correlation 
(equation 2.20) in Figure 6.16 where the latter is seen to overpredict the mean 
Nusselt number. This difference can be attributed to several factors. A small 
diameter tube with a larger diameter smooth entrance length could cause severe 
secondary flows. Also, because of the smaller diameter tube a different size twisted 
tape (twist ratio, y=3) was used, which had a slightly looser fit than the tape used 
with the smooth tube. Additionally, Prandtl number effects need to be investigated 
more closely, as Figure 6.16 shows a definite trend with inlet temperature. 

3. Cu/Ni Finned Tube - HEATEX Insert 

Data for the Cu/Ni finned tube using the HEATEX insert are presented 
in Figures 6.17 and 6.18. Figure 6.17 (Nu vs. Re) again shows a dependence of the 
data on the coolant inlet temperature. Since the HEATEX element used for this 
tube was a smaller diameter (10.16 mm) with a different loop density than that used 
for the smooth tube, a new correlation needed to be developed. Again, using the 
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form of the correlation proposed by Oliver and Aldington, the final proposed correlation 
is given by: 

Nu m - Q. 126 Pr aM Re a6S1 (6.6) 

The data and resulting correlation are compared to the correlation developed 
by Mazzone (equation 2.24) for this insert in Figure 6. 18. It can be seen that Mazzone’s 
correlation is adequate for predicting the inside heat transfer coefficients at low Reynolds 
number (Re < 1400), where both correlations merge. However, at higher Reynolds 
numbers, Mazzone’s correlation overpredicts the present experimental data. The 
differences in the two correlations are primarily due to the limited amount of data 
available to Mazzone when developing his correlation. Furthermore, the data he used 
were taken on a different experimental facility (by the manufacturers). As seen before, 
the data shown in Figure 6.18 still indicate a dependence on inlet temperature and hence 
the value of 0.46 for the Prandtl number exponent needs to be reinvestigated and 
optimized. 
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Figure 6.10 Cu/Ni externally finned tube - No insert (Nu m vs. Re) results. 
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Figure 6.11 Cu/Ni externally finned tube - No insert (Nu m vs. Re) laminar flow 
results. 
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Figure 6.12 Cu/Ni externally finned tube - No insert (Nu m vs. Re) transitional flow 
results. 
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Figure 6.13 Cu/Ni externally finned tube - No insert (Nu m vs. x*) results. 
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Figure 6.14 Cu/Ni externally finned tube - No insert laminar flow correlation results. 
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Figure 6.15 Cu/Ni externally finned tube - Twisted tape (Nu,,, vs.* Re) results. 
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Figure 6.16 Cu/Ni externally finned tube - Twisted tape correlation results. 
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Figure 6.17 Cu/Ni externally finned tube - HEATEX (Nu m vs. Re) results. 
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Figure 6.18 Cu/Ni externally finned tube - HEATEX correlation results. 



83 



D. KORODENSE TUBE 



1. Korodense Tube - No Insert 

The final tube tested was the Korodense tube. This tube was first tested 
with no insert The results for the mean Nusselt number as a function of Reynolds 
number are shown in Figure 6.19 and in more detail for laminar and transitional flow 
in Figures 6.20 and 6.21, respectively. As was seen in the smooth tube case, a 
transition to turbulent flow occurred at a Reynolds number of approximately 2300. 

Figure 6.22 compares the data to Hausen’s correlation (equation 2.10) and 
the Shah and London correlation (equation 2.14). As before, the transitional flow 
data are seen as two distinct lines for data taken at +20°C and +10°C nominal inlet 
temperatures. Under laminar flow conditions, the data are seen to fall very closely 
to Shah and Londons’ correlation, but significantly above Hausen’s correlation. This 
was also found for the other two tubes. This result was expected since a mildly 
corrugated wall is not good in enhancing laminar flow and the tube behaves 
essentially like a smooth tube. Significant enhancement for Korodense tubes is not 
expected until Re D > 10 4 [Ref. 21]. 

The resulting correlation for laminar flow data is shown on Figure 6.23 
and given by: 



Nu m = 2.05 (x*)" 0 - 526 (6.7) 

The Hausen correlation (equation 2.10) used by Mazzone again under predicts the 
present measured inside heat transfer coefficient. It was expected that the mean 



84 



Nusselt number would be the same (or maybe slightly higher) as for the smooth tube 
case for the reasons already mentioned. However, the Korodense data fall slightly 
below the smooth tube data and closer to Shah and Londons’ smooth tube 
correlation (equation 2.14). No explanation can be given for this result. Data scatter 
explanations given for the other tubes also apply for the Korodense tube. 

2. Korodense Tube • Twisted Tape Insert 

The results for the Korodense tube with a twisted tape insert are shown 
in Figures 6.24 and 6.25. The same twisted tape that was used for the smooth tube 
was used for the Korodense tube. Note that the Korodense tube nominal inside 
diameter (13.35 mm) is slightly larger than the smooth tube inside diameter (13.26 
mm). This small difference in inside diameters (together with corrugation in the 
Korodense tube) led to a looser fit for the twisted tape insert. 

The mean Nusselt number as a function of Reynolds number is shown in 
Figure 6.24 which again clearly indicates the inlet temperature dependence. The data 
were again correlated using the form of Hong and Bergles’ correlation. The present 
data were compared to their correlation (equation 2.20) in Figure 6.25. The 
correlation developed for the present data is: 

Nu m = Q.569Pr 03S (Re Jy) om (6:8) 

This correlations falls below the Hong and Bergles correlation (equation 6.20) as well 
as for the smooth tube. This result is thought to be primarily due to the looser fit of 
the twisted tape inside the Korodense tube. The corrugations also leave areas of no 
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contact between the tape and the inside wall of the tube. Stagnant flow or even air 
bubbles may form in these gaps, especially at low Reynolds numbers, causing an 
interruption of the swirled flow. The thermal boundary layer, therefore, is not as 
efficiently mixed as in the smooth tube. For these reasons, lower inside heat transfer 
coefficients than those assumed by Mazzone should be expected. 

3. Korodense Tube - HEATEX Insert 

The last tube-insert condition tested was the Korodense tube with the 
HEATEX insert The same HEATEX insert used for the smooth tube was used 
here. Figure 6.26 shows the mean Nusselt number as a function of Reynolds number 
and again shows an inlet temperature dependence. The resulting correlation for the 
present data is given by: 

Nu = 0.282 p r 0A6 Re 0A " (6.9) 

fli 

The data and this correlation are compared with the correlation Mazzone used 
(equation 2.23) in Figure 6.27. It is evident that Mazzone’s chosen correlation 
overpredicts the inside heat transfer coefficient significantly. This again is thought 
to be primarily due to the incomplete mixing of the thermal boundary layer caused 
by non-contact regions existing between the tube wall and HEATEX petals. 
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Figure 6.19 Korodense tube - No insert (Nu m vs. Re) results. 
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Figure 6.20 Korodense tube - No insert (Nu m vs. Re) laminar flow results. 
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Figure 6.21 Korodense tube - No insert (Nu,,, vs. Re) transitional flow results. 
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Figure 6.22 Korodense tube - No insert (Nu m vs. x*) results. 
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Figure 6.23 Korodense tube - No insert laminar flow correlation results. 
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Figure 6.24 Korodense tube - Twisted tape (Nu m vs. Re) results. 
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Figure 6.25 Korodense tube - Twisted tape correlation results. 
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Figure 6.26 Korodense tube - HEATEX (Nu m vs. Re) results. 
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Figure 6.27 Korodense tube - HEATEX correlation results. 
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E. INSERT COMPARISON BETWEEN TUBES 



Figures 6.28, 6.29 and 6.30 each compare the three correlations for each tube 
for no insert, twisted tape, and HEATEX inserts, respectively. It can be seen that 
with no insert (Figure 6.28), the Korodense correlation is about 10% below the 
smooth tube data. Although one may expect no increase in the heat transfer with a 
roped tube under laminar flow conditions, one certainly would not expect to find 
such a decrease with a similar diameter tube. This difference, therefore, can only be 
attributed to experimental uncertainty, probably in the wall temperature 
measurement The finned tube correlation is lower than both the Korodense (— 5 %) 
and smooth tube (=*15%) correlations; this tube had a smaller inside diameter and 
it may be that different induced flow patterns are set up causing a reduction in heat 
transfer. Experimental uncertainty must also not be ruled out as to the cause of the 
difference. 

Figure 6.29 shows the same three tube correlations, but for the twisted tape 
insert The smooth tube again has the highest performance. The finned tube is =20% 
below at a high Re number and nearly 40% lower at a low Re number. This could 
be due to a different twist ratio used for the finned tube due to the smaller inside 
diameter. The Korodense tube is similar to the finned tube, giving slightly lower heat 
transfer at high Re number and slightly higher heat transfer at low Re. The reason 
for this tube being significantly lower than the smooth tube is though to be due to 
the tape not contacting the inner surface of the tube as well as for the smooth tube. 
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thereby leaving regions of stagnant flow. These regions would tend to inhibit mixing 
and lead to a lower heat transfer performance. 

For the HEATEX comparison (Figure 6.30), similar conclusions can be made, 
although the finned tube seems to have better relative performance indicating that 
the small diameter HEATEX insert used in the small diameter tube is more effective 
than the twisted tape. The difference between the smooth and Korodense tubes are 
of similar magnitude to that shown in Figure 6.29 (about 40% at high Re). 



97 



NO INSERT COMPARISON 
UNDER LAMINAR FLOW 




Pr=70 



Figure 6.28 No insert comparison under laminar flow for all three tubes. 
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Figure 6.29 Twisted tape insert comparison under laminar flow for all three tubes. 



99 



200 



HEATEX INSERT COMPARISON 
UNDER LAMINAR FLOW 




0 500 1000 1500 2000 2500 

Reynolds 

Nu vs Re 
Pr = 70 



Figure 6.30 HEATEX insert comparisons under laminar flow for all three tubes. 
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F. LAMINAR FLOW ENHANCEMENT RATIO 



The comparisons made in Figures 6.28 through 6.30 are made more 
complicated by the fact that the uncertainty associated with the wall temperature 
measurement is unknown and most likely different for each tube. For this reason, it 
is not known how much of the divergence seen in these figures is due to differences 
in thermocouple attachment (giving uncertainty to the wall temperature 
measurement) and how much is an actual flow induced phenomenon. A comparison 
of the relative effects of the three inserts for a given tube is more accurate, however, 
due to the fact that any uncertainty for a given tube wall thermocouple attachment 
combination would be the same. 

To determine the degree of tubeside heat transfer enhancement achieved by 
the use of the twisted tape and HEATEX inserts under laminar flow, an 
enhancement ratio of the mean Nusselt number for the insert condition (Nu^,,) to 
that for the no insert condition (Nu noiniert ) was calculated. The enhancement ratio is 
therefore defined as: 



€ = 



Nu. 

insert 

Nu . 

no insert 



( 6 . 10 ) 



Figures 6.31 through 6.33 present mean Nusselt number as a function of Reynolds 
number under laminar flow conditions for each tube under the three insert 
conditions. For the purposes of this study, enhancement ratios were calculated at the 
two extreme Reynolds numbers of 200 and 2300 with a constant Prandtl number of 
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70 using the correlations developed in the previous sections. Table 6.2 summarizes the 

enhancement ratios for each tube and insert condition. 

It is evident from these figures that the greatest enhancement is achieved using the 

HEATEX insert. An enhancement factor of over seven was achieved for the smooth tube 

using the HEATEX insert at Re = 2300. This reduced to about 2.6 at low Re. 

Significantly smaller enhancements were achieved for the twisted tape inset for the 

smooth tube, indicating that for laminar flow, the stripping and mixing effect of the 

HEATEX element is considerably more effective than the simple swirl flow effect 

provided by the twisted tape. It should be mentioned here that HEATEX was developed 

specifically with laminar flow enhancement in mind. In turbulent flow (Re > 10*), the 

enhancements obtained are much smaller (around 1.5 to 2) and both twisted tape and 

% 

HEATEX elements give similar enhancement values. 



TABLE 6.2 SUMMARY OF LAMINAR FLOW ENHANCEMENT RATIOS 
(defined by equation 6.10) 





SMOOTH TUBE 


FINNED TUBE 


KORODENSE TUBE 


TWISTED 

TAPE 


HEATEX 


TWISTED 

TAPE 


HEATEX 


TWISTED 

TAPE 


HEATEX 


Low Re 
(Re = 200) 


1.79 


2.66 


1.49 


2.86 


1.92 


2.58 


High Re 
(Re = 2300) 


4.04 


7.19 


3.46 


6.71 


2.93 


4.25 
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For the finned and Korodense tubes, the HEATEX element also provides the 
highest enhancements. At high Re, the enhancements provided by the twisted tape 
and HEATEX inserts are lower than those obtained with the smooth tube for the 
reasons mentioned previously. However, at low Re, the enhancements for the three 
tubes are about the same for a given insert, indicating the "stagnant" regions 
mentioned earlier have less of an effect on the total heat transfer. 
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Figure 6.31 Insert condition comparison for laminar flow in a smooth tube. 
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Figure 6.32 Insert condition comparison for laminar flow in a Cu/Ni externally 
finned tube. 
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Figure 6.33 Insert condition comparison for laminar flow in a Korodense tube. 
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G. REPROCESSED REFRIGERANT CONDENSATION DATA 

Figure 6.34 shows the condensation data of Mazzone [Ref. 6] with the outside heat 
transfer coefficient recalculated using the present correlations. The data presented are for 
the top copper smooth tube. Similar recalculations can also be carried out for the other 
two tubes (Cu/Ni finned and Korodense tubes) as well as for condensation bundle data. 
The data of Mazzone are reprocessed using the Modified Wilson Plot Technique (see 
Mazzone [Ref. 6] for details). This technique essentially lets the leading coefficient of 
the correlation ’float’ to try to smooth out any uncertainty in the correlations being used. 
As previously stated, the chosen correlation for the no insert condition (even when 
reprocessing using the Wilson Plot Technique) gave values of inside thermal resistances 
greater than the measured overall thermal resistance, thereby yielding negative values of 
the outside heat transfer coefficient, h°. This explains why Figure 6.34 has no data for 
Mazzone’s no insert condition. 

It can be seen that by using the present developed correlations for the no insert 
condition, the calculated values of h„ are significantly improved (i.e. positive values). 
In addition, for the other two inserts, the new correlations increase h 0 by about 10% for 
the twisted tape insert and nearly 20% for the HEATEX insert. The newly calculated 
values of h 0 when using these inserts are now very close and just above the Nusselt 
prediction as may be expected for low vapor velocity film condensation. Mazzone 
commented that his data were below Nusselt by about 10% 
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and thought that this was a consequence of the length of his condenser tubes violating the 
uniform wall temperature assumption of the Nusselt theory. The fact that these data are 
so close gives confidence in the new correlations, as the insert itself should have no 
influence on the magnitude of h,. By significantly reducing the inside resistance, the 
insert simply increase the accuracy of h^. Of course, the no insert condition should also 
give similar values of h 0 . Figure 6.34 shows that they are too high, indicative of the fact 
that the correlation for the no insert condition, although, much improved, is still not ideal. 
If one recalls the scatter obtained in the present data for the no insert case, together with 
the fact that the inside resistance is still dominating, then the results shown in Figure 
6.34 are not too surprising and still very encouraging. It is recommended, therefore, that 
all future refrigerant condensation experiments carried out on the bundle facility should 
be conducted with the inserts using the newly formulated correlations. 

One final note is the fact that the present correlations were formulated 
using a constant wall heat flux, whereas for the condensation experiments, one 
should ideally use correlations developed under uniform wall temperature conditions. 
One may expect that the inside heat transfer coefficient would increase slightly under 
constant wall heat flux conditions, thereby decreasing the outside heat transfer 
coefficient. However, from Figure 6.34 it is clear that the use of the constant wall 
heat flux correlations seems to increase h 0 . Thus it appears that with the inserts, the 
inside resistance is so reduced that the boundary condition imposed has little 
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Outside heat transfer coeff./(W/m2K) 



influence. Future investigators should look into this boundary value problem in 
more detail. 
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Figure 6.34 Reprocessed refrigerant condensation data using new correlations. 
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VII. CONCLUSIONS 



The following conclusions are based upon results obtained in this study. 



A. GENERAL 



1. A test apparatus was successfully designed and built to accurately determine 
inside heat transfer coefficients for the same tubes and insert conditions used 
in the condensation experiments. 

2. Correlations to accurately determine inside heat transfer coefficients under 
laminar flow conditions both with and without inserts were successfully 
developed. 

3. Reprocessing of Mazzone’s [Ref. 6] refrigerant condensation data using these 
new correlations has yielded more accurate values of the outside heat transfer 
coefficient (Note: the name of the program used for the reprocessing is 
DRPCON7.) 



B. SPECIFIC 



1. A departure from laminar flow into a transitional flow exists at Reynolds 
numbers of approximately 2300 for the smooth and Korodense tubes and 2700 
for the Cu/Ni finned tube. This contradicts Mazzone’s [Ref. 6] assumption of 
laminar flow' up to a Reynolds number of approximately 4000. 

2. Laminar flow results for the no insert condition were in close agreement with 
the Shah and London [Ref. 9] correlation for a constant wall heat flux 
boundary condition for all three tubes tested. 

3. The Hausen [Ref. 10] correlation used by Mazzone [Ref. 6] significantly 
overpredicted the inside thermal resistance and led to nonsensical results for 

ho- 
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4. The smooth tube/twisted tape data were in close agreement with the Hong 
and Bergles [Ref. 13] correlation. Since this correlation was used by Mazzone 
[Ref. 6], his results for outside heat transfer coefficient should be fairly 
accurate for this tube. 

5. For the Cu/Ni finned and Korodense tubes with the twisted tape insert, the 
Hong and Bergles correlation predicted higher inside heat transfer coefficients 
than the experimental data over the entire Reynolds number range tested. 
This could have been due to a different twist ratio (for the Cu/Ni finned tube) 
and stagnant regions (for the Korodense tube). 

6. Experimental data for all three tubes using the HEATEX insert yielded 
• results that were well below the correlations developed and used by Mazzone 

[Ref. 6] over the entire Reynolds number range tested. 

7. HEATEX inserts provided the best inside heat transfer enhancement for all 
three tube. Enhancement ratios of over a factor of seven were found. 



Ill 



VI I L RECOMMENDATIONS 



Based on the results of this study, the following recommendations are made: 



1. Further refinement of the correlations to better account for different coolant 
inlet temperatures should be made by optimizing the Prandtl number 
exponent 

2. These correlations should be used to reprocess all refrigerant condensation 
data and should be used in future refrigerant condensation experiments 
conducted at NPS. Note that since these correlations were developed by 
omitting coolant velocities of 0.2 and 0.3 m/s, then the condensation data 
should do the same. 

3. HEATEX inserts provide the best enhancement of tubeside heat transfer. 
Due to the data reduction technique used, HEATEX inserts should therefore 
be used in all tubeside refrigerant condensation experiments to increase the 
accuracy of the outside heat transfer coefficient 

4. Transitional flow correlations can also be developed for the no insert 
condition since condensation experiments typically move into this flow region. 
However, since the transition region is not generally well developed, any 
developed correlations for this region will be useful only for this particular 
apparatus. 



112 



APPENDIX A. FLOW METER CALIBRATION 



The flowmeter used to measure flow through the instrumented tube was 
initially calibrated at 0°C and 24°C by Mazzone [Ref. 6]. However an additional 
calibration temperature of -15°C was needed for this study. The flowmeter was 
calibrated using the same procedure as Mazzone and is summarized as follows. 

The outlet of the flowmeter was disconnected from the test section and 
connected to an empty 55 gallon drum. The drum was placed on a scale with a 
resolution of 0.5 lbs. The coolant pump was started with the flow through a bypass. 
The flowmeter was initially opened to maximum flow and then throttled back to 
10%. The weight of fluid in the drum was recorded and the timer started. The weight 
was then recorded at the end of the calibration time. This procedure was repeated 
in 10% increments up to 100% flow. The calibration run was then repeated from the 
maximum flow back to 10% flow in 10% increments. The mass flow rate was then 
calculated in kg/s based on an average of the two readings. The data for the 
flowmeter calibrations as a function of nominal flow rate are shown in Table A.l. 

A least squares linear regression was performed for the flowmeter at each of 
the three temperatures for which the calibration runs were carried out The linear 
regression equation takes the following form: 
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m - a+bN 



(A.1) 

where a is the y-intercept in kg/s, b is the slope in kg/s/% and N is the flowmeter 
setting in %. The results are shown in Table A.2. 



TABLE A.1. MASS FLOW RATE (KG/S) AT -15, 0 AND 24°C 



FLOW METER A 


NOMINAL FLOW 


MASS FLOW RATE (kg/s) 


RATE 








(%) 


-15°C 


o°c 


24°C 


10 


0.0114 


0.0149 


0.0172 


20 


0.0295 


0.0363 


0.0403 


30 


0.0483 


0.0544 


0.0621 


40 


0.0679 


0.0794 


0.0851 


50 


0.0866 


0.1003 


0.1120 


60 


0.1078 


0.1179 


0.1339 


70 


0.1258 


0.1388 


0.1556 


80 


0.1453 


0.1563 


0.1769 


90 


0.1613 


0.1846 


0.2028 


100 


0.1824 


0.2019 


0.2151 



TABLE A-2. FLOWMETER CALIBRATION REGRESSION RESULTS 



FLOW METER A 




-15°C 


o°c 


24°C 


Slope 


0.0019 


0.00206 


0.00225 | 


Y-Intercept 


-0.0098 


-0.00414 


-0.00329 



In order to obtain flowrates at any intermediate temperatures, a simple 
interpolation procedure was incorporated into the data reduction program 
"DRPSING." The results of the calibrations for this flowmeter are shown in Figure 



A.1. 
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FLOWMETER 

CALIBRATION @ -15, 0, 24C 



Flow Rate (kg/s) 




Figure A.l Flowmeter Calibration. 
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APPENDIX B. HEAT LOSS TO AMBIENT EXPERIMENT 



A series of heat loss experiments were conducted to account for the amount 
of heat lost to ambient To accomplish this, the coolant was drained from the test 
tube and the inlet and outlet capped off to prevent any air flow through the tube. 
Heater power was set at a desired level and wall temperatures, monitored by the 
DAU. Average wall temperatures were determined every minute for the first 15 
minutes and then every four minutes until an equilibrium temperature was reached 
(after approximately 60 minutes). At that point, any additional heat added could be 
considered lost to ambient This procedure was repeated at five nominal power 
settings, 1W, 3W, 5W, 7W, and 9W. All data were taken at an ambient temperature 
of 22.0°C. 

Figure B.l is a plot of average wall temperature versus time for the 1W and 
7W power settings. It can be seen that an equilibrium average wall temperature was 
reached after approximately 60 minutes for these two power settings (as it was for 
the other three power settings.) This graph shows that if the average wall 
temperature in the experiments reached 46°C or 74° C, then the heat lost to the 
atmosphere was 1W or 7W, respectively. 

Table B.l shows the data for the equilibrium average wall temperature as a 
function of the actual power setting (q Upe ). The data are shown graphically in Figure 
B.2. A least squares linear regression was performed and the resulting equation was: 
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(B.l) 



_ (7V-44.15rC 
qtoa 4.39 W/°C 

Therefore, for a given average wall temperature, the amount of heat lost to 
ambient could be determined. Equation B.l was incorporated into the data reduction 
program "DRPSING." 



TABLE B.l EQUILIBRIUM AVERAGE WALL TEMPERATURE RESULTS 



NOMINAL POWER 


ACTUAL POWER (W) 


EQUILIBRIUM T w (°C) 


1W 


0.9412 


46.47 


3W 


3.122 


59.23 


5W 


5.126 


67.90 


~7W 


6.790 


74.34 


9W 


9.560 


84.85 
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HEAT LOSS TO AMBIENT 

Temperature vs Time 
1W and 7W 



Temperature (C) 




40 60 

Time (min) 



100 



Tamb ■ 22.0 

Figure B.l Heat loss to ambient; Average wall temperature versus time for 1W 
and 7W. 
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HEAT LOSS TO AMBIENT 

Temperature vs Heat Input 



Temperature (C) 




Tamb • 22.0 



Figure B.2 Heat loss to ambient; Average wall temperature versus heat input. 
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APPENDIX C. FRICTIONAL/HEAT LEAKAGE EXPERIMENT 



As coolant flows through the thermocouple wells/mixing box and along the 
tube, a small temperature rise results from frictional dissipation. When the coolant 
inlet temperature drops below ambient, an additional temperature rise can be 
attributed to heat leakage into the coolant from ambient. Since there was no means 
available to determine these terms separately, they were considered as one heat 
leakage term (q^). 

To determine this heat leakage term, the power to the heating elements was 
secured and coolant passed through the tube being tested at a given flowmeter 
setting. Any temperature rise detected was therefore due to frictional effects or heat 
leakage from ambient This temperature rise was determined at flowmeter settings 
of 15% and 100% for each tube and insert condition at all four nominal inlet 
temperatures. 

Data showed that flowmeter setting, tube and insert condition had no 
substantial effect on the temperature rise and that only the coolant inlet temperature 
made any difference. It can therefore be concluded that the heat leakage term was 
primarily due to heat addition from ambient. The results are shown in Table C. 1. 
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TABLE C.1 HEAT LEAKAGE RESULTS 



INLET TEMPERATURE 


TEMPERATURE RISE 


+20°C 


o.o°c 


+ 10°C 


0.07°C 


0°C 


0.08°C 


-10°C 


0. 18°C 



These corrections were incorporated into the data reduction program 
"DRPSING'' by directly subtracting this temperature rise from the measured outlet 
temperature. 
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APPENDIX D. POWER CALIBRATION 



Actual heating element power (and therefore heat flux) was based upon the 
applied voltage and total heating element resistance. To measure voltage and 
current, the Data Acquisition Unit made use of a voltage and current sensing unit 
However, a power factor needed to be applied to these sensed values to determine 
actual values required for the power/heat flux calculation. 

Voltage and current meters were placed in the circuit downstream of the power 
source and prior to the first heating element as shown in Figure D.l. The voltage 
rheostat was varied over its entire range and voltage and current sensed by the 
meters and in the DAU were recorded (Table D.l). For each setting, actual power 
was calculated by: 



P = VI 

and total heating element resistance by: 




(D.l) 



(D.2) 



In addition, voltage (PF V ) and current (PF,) power factors were determined by: 
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PF 



V 

ACTUAL 

~v 

CH29 



and PF, = 



l ACTUAL 
^CH25 



(D.3) 



where the subscripts "CH25" and "CH29" refer to the DAU sensing unit channel 



assignments and "actual" to the meter readings. All the above readings and 
calculations were done for five rheostat settings as shown in Table D.l. 

An average heater resistance (R), current (PF!) and voltage (PF V ) power factors 
were then calculated. Heater power could then be calculated by one of the following 
forms of equation: 



P = (PF ,)(PF v ) V cm9 



(D.4) 



or 



P * (TF.VamfR < D ' 5 > 

The current power factor varied too much to be utilized in a simple calculation. 
Therefore, it was decided to make use of equation (D.5). The final equation to 
determine the average heating element power (q Upe in Watts) for a given rheostat 
setting was: 

4*. - [59.66^/49.6 

This equation was incorporated into the data reduction program "DRPSING." 
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TABLE D.l POWER CALIBRATION RESULTS 



POWER 

(W) 


V 

(VOLTS) 


I 

(AMPS) 


R 

(OHMS) 


^CH» 

(VOLTS) 


IcHU 

(AMPS) 


PFy 




3131 


30.7 


0.61 


5032 


03135 


0.0106 


59.78 


57.50 


76.14 


61.4 


134 


4931 


1.032 


0.0430 


59.49 


28.80 


30330 


122.5 


2.48 


4939 


2.054 


0.1050 


59.63 


23.60 


75036 


192.4 


3.90 


4930 


3.226 


0.1746 


59.64 


2230 


| 113.84 


234.0 


4.76 


4931 


3.915 


0.2170 


59.77 


21.90 



POWER CALIBRATION 
ELECTRICAL SCHEMATIC 




Figure D.l Electrical schematic of power calibration experiment. 
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APPENDIX E. TUBE AND INSERT SPECIFICATIONS 



TABLE E.1 TUBE CHARACTERISTICS/SPECIFICATIONS 



PHYSICAL 

CHARACTERISTICS 


Cu 

SMOOTH 


Cu/Ni 

FINNED 


Cu/Ni 

KORODENSE 


Inside Diameter (Dj) 


13.26mm 


10.16mm 


13.35mm 


Outside Diameter (D 0 ) 


15.88mm 


14.00mm 


15.85mm 


T/C Diameter (D r/C ) t/c 


14.36mm 


12.62mm 


14.53mm 


Wall Thickness 


1.31mm 


1.22mm 


1.25 mm 


Length 


Overall Length (L) 


152.4cm 


153.7cm 


152.7cm 


Heated Length (L H ) 


140.6cm 


143.0cm 


142.4cm 


Heated Finned Length 

(L f ) 


- 


122.3cm 


- 


Smooth Section Length 

(4) 


- 


15.7cm 


- 


T/C Location (1) 


Pos 1 


15.2cm 


17.8cm 


15.6cm 


Pos 2 


55.9cm 


57.2cm 


56.2cm 


Pos 3 


96.5cm 


96.5cm 


96.8cm 


Pos 4 


137.2cm 


135.9cm 


134.5cm 


T/C X/D for L H (2) 


Pos 1 


6.77 


2.05 


7.91 


Pos 2 


37.41 


40.80 


38.35 


Pos 3 


68.06 


79.55 


68.79 


Pos 4 


98.71 


118.27 


96.99 


Inside Surface Area for L H 


A, 


0.0586m 2 


0.0340m 2 (3) 


0.0600m 2 


T/C Attachment Method 


Silver 

Solder 


Soft Solder 


Tach Weld/ 
Liquid Metal 


(1) measured from tube inlet; (2) measurec 


i from start of heated length; (3) tor 



finned heated length. 
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TABLE E.2 INSERT CHARACTERISTICS/SPECIFICATIONS 



TWISTED TAPE 


NARROW 


WIDE ! 


thickness, 8 


0.559mm 


0.559mm 


width, w 


10.16mm 


13.26mm 


180° pitch, H 


3.05cm 


5.30cm 


twist ratio, y 


3 


4 


material 


brass 


brass 


HEATEX 


SMALL 


LARGE 


diameter 


10.16mm 


13.26mm 


material 


stainless steel 


stainless steel 
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APPENDIX F. UNCERTAINTY ANALYSIS FOR THE INSIDE 
HEAT TRANSFER COEFFICIENT 



Data file SMNI03 (smooth tube, no insert, + 10°C/10kW/m 2 ),was chosen to 
conduct an uncertainty analysis of the experimental determination of the inside mean 
heat transfer coefficient, h ; . In order to assess the effects of coolant velocity on the 
uncertainty in the measurement of h ; , this uncertainty analysis was performed at 
coolant velocities of 0.4 and 1.4 m/s (for the same data set). The calculations below 
are shown for the 0.4 m/s coolant velocity. Table F.l summarizes the results of the 
uncertainty analysis for both 0.4 m/s and 1.4 m/s coolant velocities. 

The uncertainty analysis utilizes the procedure suggested by. Kline and 
McClintock [Ref. 24]. This procedure states that if 



R=R (x^ 2 s v x n ) 

then the uncertainty in R, namely, 8R is given by the equation: 



SR 



l\ 



„ \ 2 
dR s 
— Sx 



\ 2 

dR s ' 

dx 2 





' dR~ } 

At 


2 


+ 






dx n 
n 1 





where x„ is the measured variable and 5x„ is the uncertainty in the measured variable. 

The uncertainty in the coolant temperature measurements is due to uncertainty 
in the voltage measured by the thermocouples. For the thermocouples used, this 
value was assumed to be 4 microvolts or 0.1°C. The uncertainty in A T is found from 
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the uncertainty in the coolant inlet and outlet temperatures, T a and 7 and can be 
calculated as follows: 



T . 71(14) +7(15) 
" 2 



Taking the partial derivatives we obtain: 

ar dT i 

CO CO _ A 

37(14) " 37(15) " 2 

The uncertainty in the measurement of 7 „ is thus: 



dT 



1 0-5 



N 


2 


\ 


2 


dT 




dT 




CO 


+ 


CO 




37(14) J 




37(15) j 





Based on the thermocouple uncertainty, 8T = 0.1°C, we obtain: 



57 = 0.071 °C 

CO 

The coolant inlet temperature is calculated from a single thermocouple. The 
uncertainty in this measurement is calculated as: 



37(13) 



= 1 
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and 



ST. = 

CJ 



1 o J 



dT. 
a -ST 



\ 2 



^7(13) 



/ . 



= o.rc 



or. 



Since, 



ST. = arc 

Cl 



&T.\ 



[d*T) 



=1, 



'dT^ 



[dAT) 



-1 



then the uncertainty in AT is calculated from: 



sat - [(srj 2 *(srjf J 



or. 



SAT = 0.1225°C 

For a measured AT of 3.7°C, the uncertainty is therefore: 

= 0.0332 = 3.32% 

AT 
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Coolant velocity is calculated from: 




m 



PAi 



For a coolant velocity 0.4 m/s, m was calculated to be 0.05954 kg/s. The uncertainty 
in the reading of the flowmeter is defined as the scale interpolation. This value 
corresponds to one-half the value of the smallest marked increment In this case, the 
scale interpolation factor, 8 sif is 0.5%, which corresponds to 0.00101 kg/s resolution. 
In addition, a time-wise jitter, 5^ of 0.5% was noted during the runs. Therefore, the 
uncertainty in m is given by: 



The uncertainty in the cross-sectional area of the tubes was estimated from: 



where SD; is given as 0.1 mm based on tolerances supplied from the manufacturer. 




For a the measured coolant mass flow rate: 



_ = 0.0239 = 2J9% 
m 



1 
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Therefore, for an inside diameter of 13.26 mm , 



SA 



-L = 0.0107 = 1.07% 



Assuming that the uncertainty in the evaluation of the fluid properties is negligible, 
the uncertainty in the coolant velocity is now given by: 



SV 



or for a coolant velocity of 0.4 m/s, 



~\ OJ 



8m 



m 



SA, 



Am 



SV 



— * 0.0262 = 2.62% 



The uncertainty in the calculation of coolant Reynolds number is given by: 



8Re e 

~rT 



t -OJ 



sr/ \ 2 / cn \ 2 



SV 

V 
e / 



\ D ‘, j 



or substituting appropriate numbers we find 



SRe 



Re. 



2 = 0.0282 = 2.82% 



The uncertainty in the calculation of the heat transfer rate to the coolant is 
related to the uncertainty in the measurements of coolant mass flow rate and coolant 



131 



temperature rise where: 



or 




n oj 



satV 

A7 



ii = 0.0409 = 4.09% 

q 



The" uncertainty in the heat flux is then: 



T 0-3 



Sq- 






2 


q" ' 


U ) 


a ) 





where the uncertainty in A, is given by: 




SD : 



' 7 




OJ 



Here, the uncertainty in the length of the tube, SL, is assumed to be 5 mm while SD ; , 
is assumed to be 0.1 mm (again from given manufactures tolerances). Substitution 
of the appropriate numbers yields: 



8q!' 



<h 



= 0.0147 = 1.47% 
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Finally, since h; is given by: 




LMTD 



the uncertainty in the LMTD must be calculated. The LMTD is given by: 

LMTD = — 

In 

In order to calculate the uncertainty in the measurement of the LMTD, the following 
equations were used: 



A T 



T -T . 

Wt CJ 

T ,-T 

Wt CJ> 



dLMTD 

dAT 



In 



T T . 

wt a 

T -T 

wt CO 



dLMTD 



^[(T^TJ-(T^T d )] 





1 

i . 
H 




In 


wt a 




T -T 






Wt CO 





(T^TJ^-TJ 



dLMTD 



A T 





— i 

i 

1 




In 


wt a 




T -T 






Wt CO 


m 



(T -T ) 

\ wt coJ 
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dLMTD 

dT. 



AT 





i 

i 

K 
1 




In 


m a 




T .-T 




m 


wt CO 








The uncertainty in the LMTD is therefore given by: 



SLMTD = 



(blmw 8AT )\ 

dAT 



\ 2 

dLMTD 
dT * 






) 



2 

dLMTD ST 

-w- ST “ 






) 



dLMTD £ T X 
dT. 






/ 



Before the uncertainty in LMTD can be calculated, the uncertainty in the 

measured average inside wall temperature must be determined. The average inside 

wall was determined from a total of twelve thermocouples, three at each of the four 

longitudinal positions. The uncertainty in a single wall temperature measurement is 

* 

due to the uncertainty in the voltage measurement of the thermocouple (ST = 0.1°C) 
and the uncertainty due to the influence caused by the various attachment techniques 
used. Obviously, the variations seen in the wall thermocouple readings were due to 
this latter uncertainty and some account has to be made of this as this is by far the 
dominant uncertainty. If the average temperature at a given longitudinal position is 
assumed to be the required temperature, then the standard deviation of the three 
thermocouple readings at that position could be assumed to be this latter uncertainty. 
Therefore, the uncertainty at a given longitudinal position can be represented by: 
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ST,, = — S.ST 

where i = 1,2, 3, 4 for each of the four thermocouple longitudinal locations and 



iM is the standard deviation at each of these locations mentioned above. 
dT 



Since the average inside wall temperature at each position increases along the 
length of the tube, the average inside wall temperature for the entire tube must be 
thought of as the average of the four longitudinal positions’ average temperatures 
(Tpoj(,-))- The average inside wall temperature is then determined by: 

J __ ^ pasl pos2 + ^ pos3 + T pari 

4 

Taking the partial derivatives gives: 



. «•„. _ a7V _ 1 

Therefore, the uncertainty in the average inside wall temperature measurement is: 



ST = 



dT . 



dT 



-ST 



pasl 



pasl 



T 0.5 



dT \ 

w. XT 


2 

+ 


( dT . 1 


2 

+ 


(dT \ 


2 


dT . 

pos2 j 




dT ; ^ 




dT 9014 

pos4 y 
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Considering position 1, the average inside wall temperature is: 



T T(10) + r(ll)+r(12) 
3 



For the three thermocouple readings of 

T(10) = 25.99 # C, T(ll) = 25.74*C and T(12) = 26.0l°C, 
the average inside wall temperature for position 1 is 

Tpci = 25.91“ C 

The standard deviation at this location is given by: 



S.D. = bT^ = (£ = 0.123° C 



Therefore, 



bT 



■ = ^- = 4.7 x 10‘ 3 = 0.47% 



1 posl 



Repeating this for the other three locations gives: 

bT^ = 0.497° C 

bT^ = 0.634° C 

bT po34 = 0. 145° C 



136 



Therefore, the uncertainty in the measurement of average inside wall temperatures 



is: 



57V 4(0.25)(0.123)f +[(0.25)(0.497)f +[(0.25)(0.634)f +[(0.25)(0.145)ff 3 = 0.207 

After substituting the appropriate values for 5LMTD/5AT, 3UvlTD/3T wi: 
dLMTD/dT t0 and 5LMTD/3T cJ the uncertainty in the measurement of LMTD is: 



SLMTD = 1.122°C 



and 



8LMTD_ = 1-122 = 0 0465 , 4 65% 

LMTD 24.13 

The uncertainty in the overall heat transfer coefficient, U is then given by: * 



or 






/ SLMTD \ 
[ LMTD ) 



i 

7 




[(0.0417) 2 +(0.0463) 2 



1 

7 



_1 = 0.0625 = 6.25% 
Sh t . 
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—L = 0.0625 = 6.25% 

Table F.l summarizes the results of the uncertainty analysis for the smooth copper 
tube with no insert for coolant velocities of 0.4 and 1.4 m/s, respectively at a nominal 
inlet coolant temperature of + 10°C. 

The uncertainty analysis neglects any uncertainty in the physical properties of 
the coolant For the most part, this is probably a reasonable assumption. The 
exception to this may be the viscosity due to its sensitivity to temperature. Hence, the 
uncertainty reported here is likely to be conservative. 

It is curious that the relative magnitudes of the uncertainty for the low and high 
coolant velocities are nearly the same. Two different mechanisms dominate in this 
calculated uncertainty at low and high coolant velocities. In terms of the mass flow 
rate, the time-wise jitter and scale interpolation terms are the same for both coolant 
velocities. However, the actual mass flow rate is less at the lower flow rates. Hence, 
the uncertainty in the mass flow rate calculation is more significant at the lower 
coolant velocities. At higher coolant flow rates, the uncertainty in the calculation of 
LMTD becomes the dominant term. This is because the coolant temperature rise is 
smaller at high coolant velocities but the uncertainty in the temperature 
measurements remains the same. 
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TABLE F.l SUMMARY OF RESULTS FOR SAMPLE UNCERTAINTY ANALYSIS 



Variable 


V c = 0.4 m/s 


V c = 1.4 m/s 




0.98% 


0.87% 


T* 






bT 

CO 


0.51% 


0.56% 


T 

CO 






bAT 


3.32% 


10.67% 


AT 




| 


bin 


2.39% 


0.69% 


in 








1.07% 


1.07% 


A 






sv 


2.62% 


1.27% 


e 

~v~ 

c 






8Re c 


2.32% 


1.47% 


-rT- 






ii 


4.09% 


10.69% 


4 
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5T . 

H* 


0 . 57 % 


0 . 71 % 


T 

Wt 




0 


&h_ 


4 . 17 % 


10 . 72 % 


q" 






SLMTD 

LMTD 


4 . 65 % 


16 . 32 % 


*, 


6 . 25 % 


19 . 53 % 
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i l 8SS 

ta. 

1015 

1030 

1035 

1040 

1045 

1065 

1080 

1085 

1030 

1035 

1125 

1130 

1135 

1140 

1145 

1150 

1155 

1160 

1170 

1180 

1215 

1220 

1225 

1230 

1235 

1240 

1245 

1250 

1255 

1260 

1261 

1263 

1264 

1265 

1266 

1267 

1268 

1269 

1270 

1271 

1273 

1274 

127S 

1276 

1277 

1283 

1284 

1285 

1286 

1287 

1290 

1300 

1310 

1315 

132S 

1345 

13SS 



APPENDIX G. PROGRAM LISTING 






OSE 



QRPSING 

This program collects and processes single instrumented tube da 



• CREATEO: APRIL l, 1992 <S. NENORY > 

8EEP 

PRINTER IS l 

PRINT USING *4X , " ’SELECT OPTION* * * 

PRINT USING *6X , **0 TAKING OATA OR REPRQCESING PREVIOUS OATA*** 

PRINT USING -ex,**! PUR6E FILES*** 

PRINTER IS 701 
INPUT leal 1 

IF Icall-0 THEN CALL Nam 
IF Icall-l THEN CALL Purge 
ENO 

SUB Main 
COfl /Cc/ C( 7 > 

CON /Fid/ Ift 

CON /Nus/ I in, Tsat ,Qdol ,Hnua , K f ,Rho f ,Hf g,Nu f ,Oo , I tube 
CON /Will/ Ooa< 4 ) J^Oia< 4 ) ,Kma( 4 ) , Iact ( Oroot( 4 ) ,Ld< 3 ) 

CON /Uil2/ Oeita, I sat ,Nsats ,Hod,Cia< 4 ) , Alpaa< 3 ) 

OIN Emf< 17),Tp<4 ),T< 17 ) , Ho( 4 ) ,Qdp< 4 ) ,Uo< 4 > , P c ( 4 ) , X s tar< 4 ) f Nux < 4 ) , Nus s< 4 ) 
01N Xxx l < 4 ) , Xxx2< 4 ) ,Xxx3< 4 ),Xxx4< 4 ),XxxS< 4 ),XxxG< 4 ) 

OIN Num< 4 ) ,Nuag< 4 ) ,Ldd< 4 ) , Nu< 4 ) , Xs< 4 ) 

OATA 1.0,5.172,0.226,0.063.5.172 
READ Cia< * ) 

OATA 0.10086091, 25727. 94369, -767345. 8295, 78025S9S. 81 
OATA -9247486589 , 6 . 97688EI l , — 2 . 66192613, 3 . 94078EI 4 
R6A0 C(*> 

OATA 0.015875,0.014000,0.015850.0.0159,0.0 
OATA 0.013259,0.010160,0.013350,0.01345,0.0 
OATA 386.0,42.38,42.98,21.9,0.0 
OATA 0.0005588,3 

R6A0 Ooa< • ),Oia< ♦ >,Kma< •), Delta, Hod ■ Hod-H/Oi 

IF I tube-0 THEN 

OATA 6.755,37.41,58.06,98.71 

READ Ld< • ) i L6NGTH/0 IANETER RATIO FOR 4 T/C LOCS (SNOOTH) 

L-l.406 ! Heated length of smooth tube 

ENO IF 

IF I tube* 1 THEN 

OATA 2. 05, 40.80.79.5S, 118.29 

READ Ld< • > ' FIN LENGTH/O IANETER RATIO FOR 4 T/C LOCS (CU/Nl FIN) 

L-l.223 1 Heated finned length of Cu/Ni fin tube 

END IF 

IF I tube-2 THEN 

OATA 7.908,38.35,63.792,96.987 . 

R6A0 Ld< • ) ' LENGTH/O IANETER RATIO FOR 4 T/C LOCS (KOROOENSE) 

L-l.422 * Heated length of Korodenae tube 

ENO IF 

Jset-0 

Okaccpt-l 

SEEP 

INPUT -ENTER HONTH, OATE ANO TINE (NN:00:HH NNtSS )* ,Otg« 

OUTPUT 709; *TO* .Otg« 

8EEP 
Ihard- 1 

INPUT -UANT A HAROCOPY PRINTOUT < l -OEF-YES , 0-N >* , Ihard 
8EEP 
I in-l 

IF Ihard- l THEN PRINTER IS 701 
SEEP 



141 



13S0 INPUT "GIVE A NAME FOR THE NEW OATA FILE’.FileS 

1 35 5 CREATE BOAT File$,20 

1370 ASSIGN 0Fi le TO File* 

1375 BEEP 

1380 INPUT ’ENTER TUBE COOE (0-SMOOTH CU.1-FINNEO CU/NI , 2-KOROOENSE , 3-TIT 

FIN )“ , I tube 
1385 BEEP 

1384 INPUT ’0-NO INSERT, l-TUISTEO TAPE ( UIIOE > , 2-HEATEX (L), 3-HEATEX (S) 
, 4-TUISTEO TAPE (NARROU)’, Insert 

1385 IF Insert-1 THEN Hod-4 

1386 IF Insert-4 THEN Hod-3 

1388 BEEP 

1389 INPUT ’NOMINAL SUMP TEMP. 0—10 C, 1-0 C, 2-+10 C, 3-+20 C", I temp 

1390 BEEP 

1391 Egrat-54 

1392 INPUT ’ENTER EG CONCENTRATION ( UT PERCENT ) .OEFAULT - S4X’,Egrat 

1393 BEEP 

1395 • ENTER 709OtgS 

1400 OUTPUT 0Fi le ; OtgS 

1405 OUTPUT SFile: Itube.Egrat ,0dl,0d2,0d3,0d4,0d5 

1-430 PRINT USING ’10X,”FILE NAME: ” , 12A’ ; F i le* 

1435 PRINT 

1S9S lout-1 

1600 BEEP 

1605 INPUT ’WANT TO CREATE AN OUTPUT FILE? ( 1-OEF-YES , 0-N0 ) * , lout 
1610 IF lout — 1 THEN 

1615 • BEEP 

1620 INPUT ’ENTER A NAME FOR OUTPUT FILE'.FoutJ 

1625 CREATE BOAT FoutS.S 

1630 ASSIGN 0Fout TO FoutS 

1635 ENO IF 



1640 

1646 

1647 

1649 

1650 

1651 
1653 
1655 
1660 
1661 

1665 

1666 

1667 

1668 

1670 

1671 

1672 

1673 

1674 

1675 

1676 

1677 
1679 
1679 
1690 
1631 
1682 

1633 

1634 



IF I tube-0 THEN Otc-. 01436 ' 0 IAMETER TO T/C 

IF I tube- 1 THEN Otc-. 0126 ' OIAMETER TO T/C 

IF I tube-2 THEN Otc-. 01453 ■ OIAMETER TO T/C 

00- Ooa< 1 1 ube ) 

001- Ooa( I tube ) 

Oi*Oia( I tube ) 

Km*Kma< I tube > 

Ax-PI»0l’2/4 I Cross-sectional area 

IF I tube-2 THEN Aa-. 00014 i 8ASE0 ON UOLUERINE 

Ao*PI*Oo»L 

A i»P I »Oi »L 

IF I tube-2 THEN Ao-.0S»L ' 8ASE0 ON UOLUERINE 

IF I tube-2 THEN Ai-Ao/ 1.186 ' BASEO ON UOLUERINE 
Per im-PI»0i 



BURIAL (SMOOTH) 
BURIAL (CU/NI FIN) 
BURIAL (KOROOENSE) 



OESIGN OATA 



OESIGN OATA 
OESIGN OATA 



Anarea-PI ♦( Oo ‘2-Oi’2 >/4. 



Ldd( l )-L/< 4. »Oi ) 

Ldd( 2 >-L/( 2. *Oi ) 

Ldd( 3 )-3. *L/< 4. »Di ) 

Ldd( 4 )-L/0i 
IF I tube-0 THEN 

U len 1 * . 0627 ■ UNHEATEO INLET ENO OF COFPER SMOOTH TU8E ONLY 

U len2- . 0572 ' UNHEATEO OUTLET ENO OF COPPER SMOOTH TUBE ONLY 



ENO IF 

IF I tube- 1 THEN 



U len 1 » . 045 
U len2* . 06 
Ocuo- . 015875 
Ocui-. 013259 



i UNHEATEO INLET ENO OF CU/NI FIN TUBE ONLY 
• UNHEATEO OUTLET ENO OF CU/NI FIN TUBE ONLY 
I UNHEATEO OUTSIOE OIAMETER OF CU/NI FIN TU8E 
i UNHEATEO INSIOE OIAMETER OF CU/NI FIN TU8E 
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1685 

1686 

1687 

1688 

1689 

1690 

1691 

1692 

1693 

1694 

1700 

1701 

1702 

1703 

/ 9 )? 

1704 

1705 

1706 

1707 

1708 

1709 

1710 

1711 

1712 

1713 

1714 

1715 

1716 

1717 

1718 

1719 

1720 

1721 

1722 

1723 

1724 

1725 

1726 

1727 

1728 

1729 

1730 

1731 

1733 

1734 

1735 

1737 

1738 

1739 

1740 

1741 

1742 

1744 

1745 

1746 

1747 

1748 

1749 

1750 

1751 



Lcui - .112 I HEATED INLET SMOOTH 6N0 OF CU/NI FIN TUBE 

Lcuo-.097 i HEATEO OUTLET SMOOTH ENO OF CU/NI FIN TUBE 

Per lm-PI *Ocui 

Anarea-PI ♦( Ocuo*2-Ocui “2 )/4. 

ENO IF 

IF Itube-2 THEN 

Ulen 1- . 050 ' UNHEATEO INLET ENO OF CU/NI KOROOENSE TU8E ONLY 

Ulen2-.055 ' UNHEATEO OUTLET ENO OF CU/NI KOROOENSE TUBE QNL/ 

ENO IF 

Hr-50 i HEATER RESISTANCE 

! 

PRINTER IS 1 
BEEP 

INPUT "WHAT IS REQUIREO NOMINAL COOLANT VELOCITY IN m/a (MAX. - 1.4 m 

, Vega t 

OUTPUT 709; M AR AF13 AL13 VR5“ 

OUTPUT 709; "AS SA“ 

Vaum-0 

FOR J-l TO 5 
ENTER 709 ;E 
Vaum-Vaum+E 
NEXT J 

Emf ( 13 )-Vaum/5 
T ( 13 )-FNTvav< Emf ( 13>) 

Tavgl*T< 13) 

IF Tavgl<-15 OR Tavgl>24 THEN 

PRINT "COOLANT TEMP.. NOT IN RANGE -IS TO 24 C H 
GOTO 1703 
ENO IF 

Rhoegl *FNRhoeg( Tavgl ,Egrat ) 

Mdot 1 -Rhoegl • Ax ♦Vega t 
Pea 15-( Mdot 1 + . 0098 )/ .0019 
Pca0-< Mdot 1+. 004 14 )/ .00206 
Pca24-(Mdot 1>. 00329 )/ .00225 
IF Tavgl<0 THEN 

Pc(0)-Pcal5+< (Tavgl + 15 ) • ( Pca0-Pca IS ) )/ IS 
ELSE 

Pc( 0 )-Pca0+( Tavgl *( Pca24-Pca0 ) )/24 
• ENO IF 

PRINT “SET FLOWMETER REAOINGS CORRESPONDING TO:'* 

PRINT USING “GX,"-* OF METER A - " “ , 000 . 0 “ ; Pc ( 0 ) 

PRINT 

PRINT “HIT CONTINUE WHEN REAOY M * 

PAUSE 
Repeat • 1 

! PRINTER IS 701 

I do* 1 

ON KEY 0, 15 RECOVER 1734 

PRINT USING “4X," “SELECT OPTION 

PRINT USING “GX," “0-TAKE OATA “““ 

PRINT USING *SX , “ “ 1-SET Heater Tape Heat Flux ( OEFAULT ) 

PRINT USING “6X,“ “2-CHECK FLOWMETER X SETTING 
PRINT USING “ 4X , “ “ NOTE ; KEY 0 - ESCAPE"““ 

1 Ido * deaired option 

BEEP 

INPUT Ido 

1 Set default value for input 

IF Ido>2 THEN Ido-2 
1 Take data option 

IF Ido-0 THEN 1816 
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1752 

1754. 

1755 

1756 

1758 

1759 

1761 

1762 

1764 

1765 

1766 

1767 

1768 

1769 

1770 

1771 

1772 

1773 

1784 

1785 

1789 

1793 

1794 

1795 

1796 

1797 

1798 

1799 

1800 

1801 

1803 

1804 

1805 

( 12 ) 

1806 

1807 

) 

1808 

1809 

) 

1810 

1811 

) 

1812 

1813 

1814 

1315 

1816 

1813 

1819 

1820 

1821 

1822 

1823 

1824 

1825 

1826 

1827 

1828 

1829 



IF Ido-2 THEN 1703 

Loop to set heat flu* and check temperatures 
IF Ido-1 THEN 

INPUT -ENTER DESIRED HEAT FLUX",Qdpuan 

PRINT USING " 4X t " " Qdp (wanted) Qdp (actual) Tin Delt“ , * ,, 
PRINT USING “ 4X , " " <U/pT2) (U/m*2) (C) <C> 

Read thermocouple voltages 
OUTPUT 709; "AR AF1 AL15 UR5 " 

FOR 1-1 TO 15 

OUTPUT 709; "AS SA” 

Usum-0 

FOR J-l TO 20 
ENTER 709; E 
Usum-Usum+E 
NEXT J 

Etnf ( I )-Usum/20 
T ( I )-FNTvsv( Emf ( I ) ) 

NEXT I 

Compute average temperature of coolant outlet 
Tout-( T ( 14 ) + T < 15 ) >/2 
Del t*Tout-T ( 13) 

OUTPUT 709; -AR AF25 AL2S UR5 * 

OUTPUT 709; "AS SA" 

ENTER 709 ; Amp 

OUTPUT 709; "AR AF29 AL29 UR5" 

OUTPUT 709; “AS SA" 

ENTER. 709 ; Uo 1 t 
Qt ape-( 60*Uoi t )*2/Hr 
Qdp tap e-Q tape/ Ao 
PRINT 

PRINT USING -4X,5(rtZ.3DE,3X)-.Qdpuan,Qdptape,T(13),0elt,Qtape 
PRINT 

PRINT USING " 4X , " " TEMPS AT POSITION 1 - " ,3< ODD .00 , 3X >" ; T< 10 > , T< 1 1 ) , T 



PRINT 

PRINT USING - 4X , - - TEMPS AT POSITION 2 * " , 3( ODD . 00 , 3X ) " ; T( 7 ) , T( 8 ) , T( 9 
PRINT 

PRINT USING " 4X , " " TEMPS AT POSITION 3 " - , 3( 000 . DO , 3X ) " . T( 4 ) , T ( 5 ) , T( 6 



PRINT 

PRINT USING " 4X , - - TEMPS AT POSITION 4 " - , 3( ODD . 00 , 3X > - . T< 1 > , T< 2 > , T ( 3 

WAIT 2 
GOTO 1762 
END IF 

TAKE OATA IF Im-0 LOOP 
T wm*0 
Twi -0 

PRINTER IS 701 

OUTPUT 708; -AR AF 1 AL15 UR5 " 

FDR 1-1 TO 15 

OUTPUT 709; "AS SA" 

Usum-0 

FOR J-l TO 20 
ENTER 709. E 
Usum-Usum+E 
NEXT J 

Emf( I )-Usum/20 
NEXT I 
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1840 

1845 

1850 

1860 

1865 

1870 

1871 

1872 

1873 

1874 

1876 

1880 

1881 

1883 

1884 

1885 

1886 

1887 

1888 

1880 

1891 

1893 

1894 

1895 

1896 

1897 

1898 

1923 

1925 

1926 

1930 

1935 

1940 

1941 

1942 

1943 

1945 

1946 

1948 

1949 

1951 

1952 

1953 

1954 

1955 

1956 

1960 

1970 

1990 

1996 

2000 

2005 

2010 

2015 

2020 

2025 

2030 

2035 

2040 

2045 



DATA ANALYSIS 



FOR 1-1 TO 15 

T ( I )-FNTvsv<Emf ( I ) ) 
NEXT I 

FOR 1-1 TO 12 
Twi-Twi+T( I ) 

NEXT I 



Twiav-Twi / 12 . 

T out *< T ( 14 ) + T ( IS ) )/2 
Oeit-Tout-T< 13) 



Hlc-0. ! Heat leakage correction 

IF I tube-0 THEN 

IF Itemp-0 THEN Hie-. 18 ! HEAT LEAKAGE CORRECTION AT -10 C 

IF Itemp-1 THEN Hlc-,08 ' HEAT LEAKAGE CORRECTION AT 0 C 

IF Itemp-2 THEN Hlc».07 ■ HEAT LEAKAGE CORRECTION AT +10 

IF Itemp-3 THEN Hlc-0. ! NO HEAT LEAKAGE CORRECTION AT + 

Delt-Del t-Hlc 
END IF 

Tout-Tout-Hlc 

OUTPUT 709; - AS SA” 

ENTER 709 ; Amp 

OUTPUT 709; *AR AF29 AL29 UR5" 

OUTPUT 709; -AS SA* 

ENTER 709 ; Uolt 



Qtape-( 60 # Uol t ) A 2/Hr 
Qdptape-Qtape/Ao 
Jse t-Jse t + 1 



IF Okaccpt-0 THEN Jset-Jset-1 
PRINT 

PRINT USING " 10X .""Data set number - ■*,OD*;Jset 
PRINT 

IF I tube-0 THEN PRINT USING " 10X , " * Smooth tube - " - 

IF I tube-2 THEN PRINT USING * 10X , * ‘Korodense tube 

IF I tube-1 THEN PRINT USING " 10X , - - Copper/Nickel fin-*" 

IF I temp-0 THEN PRINT USING “ 10X , " "Norn ina 1 sump temp. - -10 

IF Itemp-1 THEN PRINT USING - 10X “Nomina 1 sump temp. - 0 C -,,M 

IF I temp-2 THEN PRINT USING " 10X t “ " Nom i na 1 sump temp. - +10 C"' 

IF Itemp-3 THEN PRINT USING " 10X , ’‘“Nomina 1 sump temp. - +20 C 

IF Insert-0 THEN PRINT USING "10X,"-No insert used 

IF Insert-1 THEN PRINT USING " 10X , " "Wide twisted tape insert used ,,,, “ 
IF Insert-4 THEN PRINT USING " 1 0X f " “Narrow twisted tape insert used" 

IF Insert-2 THEN PRINT USING " 10X , “ - Large HEATEX insert used 

IF Insert-3 THEN PRINT USING ’ 10X , M - Sma 1 1 HEATEX insert used""** 

PRINT USING “ 10X "Nomina 1 set heat flux - " ** . 30E " ; Qcp tape 

PRINT USING " 10X “Nominal set coolant velocity - " *' , 0 . D “ . Uegs t 
PRINT 



Tavg-T( 13 ) +0e 1 1 * . 5 
Rhoeg-FNRhoeg( Tavg,Egrat ) 
Nuueg-FNNueg( Tavg,Egrat ) 
Mueg-Nuueg*Rhoeg 
Cpeg-FNCpeg< Tavg.Egrat ) 

Keg-FNKeg( Tavg.Egrat ) 
Preg-Cpeg*Mueg/Keg 
Mdot-FNFmcal< T( 13 ) ,Pc<0 ) ) 

Ueg-Mdot / ( Rhoeg* Ax ) 

Reeg-Ueg *0 1 /Nuueg 

Res-4*f1dot/<PI*f1ueg + (0i-4*0elta ) ) 
Qdot-Mdot*Cpeg*0elt 
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204.8 » ITERATE TO FIND INSIOE WALL TEMPERATURES 
2049 Tnl*Tujiav 

Z0S0 Tn2-Tujiav 

2051 Tnl *( Tnl + Trc2 )/2 . 

2052 Ct l-Qdo t *LOG( Otc/Oi )/< 2 . *PI *L*Km ) 

2053 Tn2*Tw iav-C t i 

2056 IF ABS( Tm2-Tnl )> . 001 THEN GOTO 2051 

2057 FOR I-I TO 12 

2058 T ( I )»T ( I )-Ct i 

2059 Twpi-Tujr+T ( I ) 

2060 NEXT I 

2061 Twn-Tujn/12 

2062 Nuuegu*FNNueg< Tun , Egrat ) 

2063 Muegw*Nuuegu*Rhoeg 

2064 TujrI ■< T< 12 )+T ( II > + T< 10 ) )/3 

2065 TuiR2-( T < 9 ) + T < 8 > + T( 7 ) )/3 

2066 T ur3*( T < 6 ) + T(S >+T( 4 ) )/3 

2067 Tuin4-(T(3)+T(2)+T(l ) )/3 

2068 TuiRi2-(T( 12 )+T( 11 )+T( 10 )+T< 9 >+T< 8 )+T< 7 ) )/6 

2069 Tujr123-<T< 12)+T<I1 ) + T( 10 ) + T( 9 > + T( 8 ) + T< 7 )+T< 6 ) + T( 5 )+T( 4 ) )/9 

2070 T1r-( T out-T ( 13 ) >/L0G< ( Tujr-T( 13 ) )/( Turn-Tout ) ) 

2071 Q1os3-(Tur- 44. l5)/4.39 ! HEAT LOSS TO ATMOSPHERE 

2072 IF Tun<45 THEN Qlo33-0 

2073 Qtape-Qtape-Qlosa 

2074 Qdpi-Qdot/Ai 

2075 Alfai-Qdpi/Tln 

2076 » FIN CORRECTION TO ACCOUNT FOR CONOUCTION TO ENOS OF TUBE 

2077 Ff i*( Kn*PeriR*Anarea ) * . 5 

2078 Alfail-Alfai 

2079 AifaxS-Alfail 

2080 Alanl*< <Alfai2*PeriR)/( Kn*Anarea ) )* .5 

2081 Ff3»Ff 1*A1 fai2" .5*( Tun-T( 13 ) ) 

2082 Ff4-Ffl*Alfai2*.S*< Tun-Tout ) 

2083 Func6»Ff 3 # FNTanh( Alanl *Ulenl ) 

2084 Func7*Ff4*FNTanh< Alanl*Ulen2 > 

2085 Func8*Al f ai 2 * A i * ( Tum-Tavg ) 

2086 Func**Func5 +Func7+Func8-Qdo t 

2087 0func6-( .S*F u nc6/Al fai2 )+2. *Ff3*Uleni/( 1 . +FNCosh( 2 . *Alanl *Ulenl ) ) 

2088 0func7*< ,S*Func7/Al f ai2 > + 2 . *F f 4*Ulen2/< 1 . +FNCosh( 2 . *Alanl *Ulen2 ) ) 

2089 0func8*Ai •( Tum-Tavg ) 

2090 0funcx*0 f uncS^Q f unc7+0f unc8 

2091 Alfail»Al f ai2-Func.</0func:< 

2092 IF A8S<AlfaiI-Alfai2>'.0S THEN GOTO 2079 

2093 • A1 f ai2*Al f a 1 1 

2095 C f eg*( Mueg/Muegu ) " . 1 4 

2096 i FINO TEMPERATURE INCREASE OUE TO UNHEATEO LENGTHS 

2097 Oe 1 t in*Func6/ < Mdot *Cpeg ) 

2098 Oel tout*Func7/ ( Mdot *Cpeg ) 

2099 Thl i n*T ( 13 KOeltin 

2100 Thl out *Tout -Oe 1 tout 

2101 Oe 1 h 1 *( Thiou t -Th 1 i n ) 

2105 IF I tube-1 THEN 

2106 Heat m-Lcui *Oelhl/(L*Lcui-**Lcuo ) > HEATEO SMOOTH INLET SECTION FOR CU 
/N I FIN TUBE 

2107 Heatout-Lcuo*Oelhl/(L+Lcui+Lcuo ) » HEATEO SMOOTH OUTLET SECTION FOR 
CU/ N I FIN TUBE 

2109 Th 1 in* Th 1 in+Heat i n 

2110 Thlout*Thlout-Heatout 

2111 Qelhl*< Thlout-Thl in ) 

2114 ENO IF 
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2115 

2116 

2117 

2118 

2119 

2120 

2121 

2122 

2123 

2124 

2126 

2127 

2128 

2129 

2130 

2131 

2132 

2133 

2134 

2135 

2136 

2137 

2138 

2139 

2140 

2141 

2143 

— Z144 

2145 

2146 

2147 

2148 

2149 

2150 

2151 

2 1S2 

2153 

2154 

2155 

2156 

21S7 

2158 

2159 

2161 

2162 

2163 

2164 

2165 

2166 

2167 

2168 

2169 

2170 

2171 

2172 

2173 

2175 

2176 

2177 

2178 



Qhl»Mdot *Cpeg*Qeihl 
Qdphl-Qhl/Ai 
Qdpxx*Func8/Ai 
Tcoi*Thlin+(Qelhl/4 ) 

Tco2-Thlin+(Qelhl/2> 

Tco3*Thlin+(3*0elhl/4) 

Ttc 1*T( 13 ) + <Ld( 0 >*0i*0elhl/l ) 

Ttc2-T< 13) + <Ld( I ) # 0i *Oelhl/L ) 

Ttc3«T( 13 >+<Ld(2 >«0i*0elhl/L ) 

Ttc4-T( 13 >+(Ld( 3)*0i*0elhl/L > 

! ALL < 1 , 2 , 3 , 4 ) 

* SECTION 1 
» SECTION 
« SECTION 

• SECTION 4 
' SECTION 1,2 

)» SECTION 1,2, 



Kegl *FNKeg< (Tcol+Thlin)/2, Egrat ) 

Keg2»PNKeg< ( Tco2+Tcol >/ 2, Egrat ) 

Keg3*FNKeg< CTco3+Tco2 )/2, Egrat ) 

Keg4-FNKeg( ( Thlout +Tco3 >/2 , Egrat ) 

Kegl2*FNKeg< (Tco2 + Thlin)/2, Egrat ) 

Kegl23*FNKeg( (Tco3+Thlin)/2, Egrat ) 

Numa ll-Alpali*Oi/Keg 
Nul-Alpl«Oi/Kegi 
Nu2*Alp2 *0i/Keg2 
Nu3-Aip3«Qi/Keg3 
Nu4«A lp4 *01 /Keg4 
Nulac 1 -A1 loc 1 *0 l/Kegl 
Nuloc2-Alloc2*0i/Keg2 
Nuioc3-Al 1 qc3*0 i /K eg3 
Nulac 4- A 1 loc 4*0 i /Keg4 
Nul2-Alpl2*0i/Kegl2 
Nul23-Alpl23«0i/Kegl23 

LOCAL NUSSELT FROM SHAH/L0N00N PAGE 128, EQN 240-242 
CALCULATE XST'AR FIRST FOR EACH THERMOCOUPLE LOCATION 
DEFINE COOLANT PECLET NUMBER 
P-eg-Pr »g*Recg 

Xs( 1 )-Ldd< 1 >/Peeg ' FOR MEAN VALUES OF h 

X s( 2 )*Ldd( 2 )/Pee g 

Xs(3)-Ldd(3)/Peeg 

X s ( 4 )*Ldd< 4 )/Peeg 

FOR I-l TO 4 

Xstarf I >-Ld< 1-1 )/Peeg ' FOR LOCAL VALUES AT T/C POSITIONS 
Nua( I )-l .302 *X star < I r< -1/3)-. 5 

DEVELOPING HYORO ANO THERMAL 80UN0ARY LAvER (S/L PU8, EQN. 253 > 
X x 1 ■ ( 1 + ( ( 220/ PI )*Xstar( I ) ) A ( -10/9 ) r < 3/5 > 

Xx2-( l + < Preg/. 0207 r< 2/3 ) )“ .5 
Xx3-PI/< 1 15.2 *Xs tar( I ) ) 



Tlmall*(Thlout-Thlin )/L06< ( Twm-Thlin )/( Twm-Thlout ) ) 
Tlml*(Tcol-Thlin )/LOG( < Twml-Thl in)/(Twml-Tcol )) 
Tlm2-( Tco2-Tcol )/LOG( ( Twm2-Tcol >/< Tum2-Tco2 ) ) 

Tlm3*( Tco3-Tco2 )/LOG< ( Tum3-Tco2 >/( Twm3-Tco3 > > 
Tlm4*(Thlout-Tco3 )/LQG( ( Tum4-Tco3 )/( Twm4 -Thlout > ) 

T 1 m 1 2 * < Tco2-Thl m )/LOG( (Twml2-Thlin)/(Twml2-Tco2 ) ) 

Tlml23-(Tco3-Thlin )/LOG< ( T w ml23-Thl in )/ ( Twml 23-Tco3 ) 

Alpall*Qdphl/Tlmall 

Alpl-Qdphl /Tlml 

Alp2-Qdphl/Tln2 

Alp3*Qdphl/Tlm3 

Alp4-Qdphl/Tln4 

Allocl*Qdphl/(Twml-Ttcl ) 

Alloc2«Qdphl/< Twm2-Ttc2 > 

A1 loc3-Qdphl/( Twm3-T tc3 ) 

Alloc 4*Qdphl/< Twm4-Tt c4 ) 

Alpl2*Qdphl/Tlml2 

Alpl23-Qdphl/Tlnl23 
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rsi in 



2179 X x. 4 * ( !+<Xx3/< Xx2*Xxl ) r <5/3 ))* <3/10 > 

2180 XxS-< l + < (220/ PI >«Xstar< I ) r<-10/9 ) r<3/10 ) 

2181 Nuss< I )-< X*4*5.364*Xx5 )-l 

2188 IF Xstar( I X .00005 THEN 

2189 Nux< I )-l ,302*Xstar< I r< -1/3 )-l 



2190 


ENO IF 






2194 


IF Xstar( I ».00l5 THEN 






2195 


Nux ( I )-4. 364 + 8. 68* < 1000 


•Xstar( I ) >* 


( - .506 ) # EXP( -41*Xstar< I ) ) 


2196 


ENO IF 






2200 


NEXT I 


t 




2202 


IF Insert-0 THEN 






2203 


PRINT USING “3X , - " 




LOCAL NUSSELT VALUES" " " 


2205 


PRINT USING "3X , " " Exp. 


Nux ( TLM ) 


Exp. Nux(8ULK) Pred. Nux< 1 ) 



Pre d. Nux( 2 ) Xstar< tc) mmm 



2206 

2207 
star* 1 ) 

2208 
star* 2 ) 

2209 

st ar( 3 ) 

2210 

st ar( 4 ) 



PRINT 

PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 



3X,--P0S 1 
3X , " "POS 2 
3X.--POS 3 
3X # * ‘POS 4 



•* ,5<MZ.3DE,5X )• ;Nul .Nulocl ,Nux< 1 ),Nuss( 1 >,X 
•• ,5<MZ.30E,SX )-;Nu2,Nuloc2,Nux<2 >,Nuss<2),X 
* - ,S( MZ.30E.5X )" ;Nu3,Nuloc3 ,Nux< 3 ),Nuss<3 ),X 
" " ,5<MZ.30E,SX )• ;Nu4,Nuloc4,Nux< 4 ) ,Nuss< 4 ) f X 



2211 

2215 

2216 
Xstar 

2217 

2218 

s t ar ( 1 ) 

2219 

st ar( 2 ) 

2220 

3 tar( 3 > 
222 1 
3t ar ( 4 ) 

2223 

2230 

2232 

2233 

2234 

2235 

2236 

2237 

2239 

2240 

2241 

2242 

2243 

2244 

2245 

2246 

2247 

2248 
22491 



PRINT 

ELSE 

PRINT USING 
at tc po s • " " " 
PRINT 

PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 



" 3X , " "POSITION t Exp. Nux ( TLM ) Exp. Nux ( 8ULK ) 



3X, ""POSITION 


1 


(LOCAL) - ’ 


3X , * "POS ITION 


2 


(LOCAL) - " 


3X, ""POSITION 


3 


(LOCAL) - * 


3X, --POSITION 


4 


(LOCAL) - " 



- ,3<MZ.30&,SX )" ;Nul .Nulocl ,X 
- , 3< MZ . 30E , 5X ) " ; Nu2 , Nu 1 oc2 , X 
- ,3<MZ.30E,5X >" ;Nu3,Nuloc3 f X 
" , 3 ( MZ . 30E , 5X ) " ; Nu4 , Nu 1 oc 4 , X 



END IF 
PRINT 

FOR 1-1 TO 4 
IF Insert-0 THEN 

Bb 1* . 0668 *Reeg*Preg/Ldd< I ) 

Bb2- 1 ♦ . 04 * < Reeg*Preg/Ldd< I ) r .6666 . 

Nueg( I )-3.66+(8bl/Bb2 ) 

ENO IF 

IF Insert-1 OR Insert-4 THEN 

Nueg* I )-Cia( Insert ) * ( l>5.484E-3 # Preg / '.7*<Res/Hcd)"1.25)*.5 
. Nuoth*. 383 *Pr eg'* .35 ♦(Res/ Ho d) A . 622 
ENO IF 

IF Insert-2 THEN 

Nueg( I )-Cia< Insert )•< Reeg" .65 )*Preg*' . 46 
ENO IF 

IF Insert-5 THEN 

Nueg( I )-Cia( Insert ) ♦( Reeg" . 76 ) *Preg~ . 46 
ENO IF 
BEEP 

CALCULATED VALUE OF INSIDE HT. TRANSFER COEFF. FROM CORRELATIONS 



2250 H i -Nueg( 4 ) *K eg/0 i 

2251' MEAN NUSSELT NUMBER FROM SHAH/LONOON PAGE 128, EQN 245-246 

2254 IF X s( I )< . 03 THEN 

2255 Num< I ) — 1 . 9S3*Xs< I ) * < - 1/3 ) 
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2056 


ELSE 








2258 


Nun( I >-4.3S4+< . 0722/Xs ( I ) ) 




2259 


ENO IF 








2262 


NEXT I 








2263 


PRINT USING 


"3X, ""POSITION 


* “ Exp. 


mean Nu Pp*d. mean Nu 


Xs tar 


at sec t ion " " “ 








2264 


PRINT 








2266 


PRINT USING 


"3X, "“POSITION 


1 (MEAN) - 


,3( MZ . 30E ,5/ )" ; Nul ,Nun( 1 ) 


, X 3 ( 1 ) 
2267 


PRINT USING 


"3X , " "POSITION 


12 (MEAN) - 


, 3< MZ . 30E ,SX )’ :Nul2,Nun( 2 


> , Xs( 2 > 
2268 


PRINT USING 


" 3X , " "POSITION 


123 (MEAN) - •* 


, 3( MZ . 30E , SX ) “ ; Nul23 ,Num( 


3 ) , Xa( 3 ) 
2269 


PRINT USING 


" 3X , " "POS ITION 


1234 (MEAN) * 


, 3( MZ . 30E.SX ) * ; Nuria l 1 , Num 


( 4 > , X s ( 4 ) 
2270 


PRINT 








2271 


PRINT 








2272 


PRINT USING "10X,""Mass 


flow rate ( kg/s ) 


- •* .MZ.30E' ;Mdot 


22731 


PRINT USING ‘10X, ""Inside Tube Die. <m> 


" “ " , MZ . 3DE“ • Oia( 1 1 ube ) 


2274 


IF Insert 


-1 THEN PRINT USING *I0X,-*t80 


OEG OVER Oia . ( H00 ) - H " , 



MZ . 30E" ; Hod 



2275 


PRINT USING 


"10X 


2276 


IF Insert*! 


OR I 


2277 


PRINT USING 


" 10X 


2278 


ENO IF 




2279 


PRINT USING 


" 10X 


2280 


PRINT USING 


" 10X 


2281 


PRINT USING 


"10X 


2282 


PRINT USING 


- I0X 


2283 


PRINT USING 


- 10X 


2284 


IF Inser t*l 


OR Ii 


18 S 


.MZ.30E" ;Res 




2285 


PRINT USING 


" i 0X 


2286 


PRINT USING 


-10X 


4 

2287 


IF lnsert*0 


THEN 


" , MZ . 30E 


" ;Nueg( 4 > 




2288 


IF Insert*! 


OR* Ii 


pe (BerglesZ )”", MZ.30E" iNuoth 


2289 


IF Insert*! 


OR Ii 


pe ( 8ergies3 > " " , MZ . 30E" ; Nueg( 4 ) 


2290 


IF Insert*2 


THEN 


.MZ.30E’ .Nuegf 4 > 




2291 


IF Inser t *3 


THEN 


"* .MZ.3DE" . Nueg( 4 > 




MMQM 

-49*. 


IF lnsert*0 


THEN 


MZ.3DE’ :Nun<4) 




2293 


PRINT USING 


'* 10X 


2294 


PRINT USING 


- 10X 


2295 


PRINT 




2296 


PRINT USING 


-3X, 


)- .T( 12 > 


,T( 11 ) , T ( 10) 




2297 


PRINT USING 


"3X, 


)- ;T(9), 


T( 8 ) , T< 7 ) 




2298 


PRINT USING 


-3X, 


>".T(S >, 


T ( 5 ) , T( 4 ) 




2299 


PRINT USING 


-3X, 


>" : T( 3), 


T( 2 >,T( l ) 




2300 


PRINT USING 


-3X, 


2301 


PRINT 





- --,flZ.30E-;T(13) 

- "" .MZ.30E" ;0elta 

* " “ ,flZ . 30E’’ : Oel t 

- ,MZ.30E";Qdot 

:U> - , MZ - 30E” ; Qtape 



' ,MZ . 3DE" . Preg 
’ , MZ. 3DE” ; Reeg 



,MZ . 30E " ; A 1 f a i 



( pos . 2 ) ( C ) 
(pos. 3) (C) 
( pos . 4 ) ( C ) 



1 , 3(000. 0D.3X 
’ ,3(000.00, 3X 
' .3(000. 00, 3X 
' ,3(000. 00, 3X 
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2302 


Nu( 1 )*Nul 






2303 


Nu( 2 )»Nul2 






2304 


Nu< 3 )*Nul23 






2305 


Nu< 4 )»Numa 1 1 






2307 


IF Insert*© THEN 






2308 


FOR 1-1 TO 4 






2310 


Xxxl( I )»LOG< Nu< I ) ) 






2311 


Xxx2< I >-LOG< Xs< I ) ) 






2312 


Xxx3< I >-LOG<Nueg< I ) > 






2313 


Xxx4< I )-LOG(Nun< I ) ) 






2314 


NEXT I 






2315 


PRINT USING "3X , “ "MEAN UALUES AT 


POS 1 


POS 1 


2 


POS 123 POS 1234" ““ 






2317 


PRINT USING “SX , * *EXP . LOG< Na ) 


,4<nZ.30E, 


4X )• iXxxl( 1 ),X 


xxi<2 ) 


, Xxx 1 < 3 ) , Xxx 1< 4 ) 






2318 


PRINT USING "SX , * “COR . LOG<Nu) <H) 


-• t 4(nZ.30E, 


4X )" ;X/.x3< 1 ),X 


xx 3(2 ) i Xaa3(3),Xxx3<4) 






2319 


PRINT USING “SX , “ “COR . LOG<Nu> < S/L ) » 


■■ ,4<nZ.30E, 


4X )" ;Xxx4< 1 ) t X 


xx4<2 ),Xxx4<3),Xxx4(4) 






2320 


PRINT USING - SX | --LOG<Xstar ) 


" " ,4<MZ . 30E t 


4X )• ;Xxx2< 1 ),X 


xx2< 2 ) ,Xxx2<3 ) ,Xxx2< 4 ) 






2321 


ENO IF 






2323 


IF Insert»2 OR Insert*3 THEN 






2324 


FOR 1-1 TO 4 






2325 


Xxxl< I )-LOG<Nu( I )/Preg" .46 ) 






232S 


Xxx2< I )”LOG( Reeg > 




• 


2327 


Xxx3( I )-LOG(Nueg( I )/Preg".4S ) 






2328 


NEXT I 






2328 


PRINT USING “3X , " “MEAN UALUES AT 


POS 1 


POS 1 


-» 

4 . 


POS 123 POS 1234"“" 






2331 


PRINT USING “6X , “ “EXP . LOG< Nu/Pr 46 ) - 


-“,4(MZ.30E 1 4X) , ‘;Xxxl(l) 1 


Xxx 1 < 2 


),Xxxl<3),Xxxl<4) 






2332 


PRINT USING “SX , “ “COR. LOG< Nu/Pr *. 4S ) * 


,4(MZ.30E 


.4X7". Xxx3< 1 ) , 


X x x 3 < 2 


) ,Xxx3< 3 ) ,Xxx3( 4 ) 






2334 


PRINT USING "SX , * “LOG< Re ) 


“" , 4( MZ . 30E 


,4X )";Xxx2< 1 ). 


Xxx2<2 


) ,Xxx2< 3 ) , Xxx2< 4 ) 






2335 


ENO IF 






233S 


IF Insert*l OR Insert»4 THEN 






2337 


FOR 1-1 TO 4 






2338 


X •. .< 1 < I )*LOG< Nu< I )/Preg" .35 ) 






2338 


Xx <2( I )*LOG( Res/Hod ) 






2340 


Xax 3< I )*LOG< Nuoth/Preg" . 35 > 






2341 


Xxx4< I ) “LOG < Nu< I ) ) 






2342- 


X x x 5 < I >»LOG< Preg*< Res/Hod )“ 1.78) 






2343 


X x x 6 < I )-LOG< Nueg< I ) ) 







2344 next I 

2345 PRINT USING 

2 POS 123 POS 

2346 PRINT USING 
Xxx 1 ( 2 ) , Xxx 1 < 3 ) ,Xxxl<4 > 

2347 PRINT USING 
Xxx3< 2 ) ,X.«.x3< 3 >,Xax3< 4 ) 

2348 PRINT USING 
X.<x2( 2 ),Xxx2(3),Xxx2(4) 

2349 PRINT 

2350 PRINT USING 
X x x 4 ( 2 ) , X x 4 ( 3 ) , X -,<4( 4 ) 

2351 PRINT USING 
X xxS < 2 ) ,Xxa 6< 3 ) ,Xxx6< 4 ) 

2352 PRINT USING 
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*3X , “ "MEAN UALUES .AT POS 1 POS 1 

1234"“" 

“6X , " "EXP . LOG< Nu/Pr- .35 ) =■ " " , 4( Ml . 30E , 4X ) " ; X/.. 1 ( 1 > , 
"SX , " “COR . LOG< Nu/Pr " .35 ) - " “ , 4( (1Z . 30E , 4X ) " , Xxx3< 1 ) , 
"6X , " “LOG< Res/y ) =■ " * , 4< MZ . 30E , 4X )" ; Xxx2< 1 > , 

"SX , " "EXP . LOG(Nu) - "“ ,4(NZ.30E,4X )" ,Xax 4( l 7, 

"SX , “ "COR. LOG(Nu) - ,4<nZ.30E,4X )" ;Xa.-. 6( l ), 

“6X,""L0G<Pr(Res/y >*1.78) * “ " , 4( MZ . 30E ,4X7", Xxx 5(1), 



,4(MZ .3DE.4X )" ;Nu 



Xx*5< 2 >,Xx*5(3),X.<x5(4) 

2353 ENO IF 

2354 IF Insert *0 ANO Reeg>2300 THEN 

2355 Nutu-.027*Reeg".8*Preg~.3333*Cfeg 

2355 PRINT USING "GX , " " TURBULENT S/T CORRELATION - 

tu 

2357 X*x7-L0G<Nutu ) 

2358 PRINT USING '*6X f " " TURB . COR. LOG<Nu) - H , 4< MZ . 3DE , 4X ) ” ; Xax7 

2359 ENO IF 

2350 Okaccpt-l 

2351 BEEP 

2352 INPUT -OK TO ACCEPT THIS SET < 1-OEFAULT-YES , 0-NO >7 - , Okaccp t 

2353 IF Qkaccpt -l THEN OUTPUT ©Fi ie s Pc< O , Enf ( * ) , Tp( * ) 

2354 IF Id-0 THEN 

2355 Okrpt-l 

2355 8EEP 

2357 INPUT "WILL THERE 8E ANOTHER RUN ( i-YES-OEFAULT , 0-NO ) - , Okrpi 

2358 IF Okrpt-l THEN 1701 

2359 ELSE 

2370 IF JseUNaets THEN 1701 

2371 ENO IF 

2372 ASSIGN ©File TO •- 

2373 IF Iout-l THEN ASSIGN ©File TO * 

2374 SUBENO 

2375 OEF FNGrad(T) 

2375 Grad— 3; 877857E-S-2*4. 71 428S7E-8*T 

2377 RETURN Grad 

2378 FNENO 

2379 OEF FNKcu(T) 

2380 ' OFHC COPPER 250 TO 300 K 

2381 Tk-T + 273 . 15 

2382 K-434-. 1 1 2 * Tk 

2383 RETURN K 

2384 FNENO 

2385 OEF FNNueg( Tc , Egr ) 

2385 • RANGE OF UALIOITY: -20 TO 20 OEG C 

2387 Tk -Tc+273 . 15 

2388 Nu I -7 . II 96507E-3-T k * ( 7 . 4863347E-5-Tk *<2. 629494 3E-7-Tk*3. 08 3332 9E -10) ) 

2390 Nu2-4 . 9237G38E-3-Tk *(4.92 I39I2E-5-Tk • ( 1 . 6437534E-7-TK ♦ l .333333 IE- 10 ) ) 

2395 Nu3-8 . 6586293E-3-Tk*(8.8837902E-5-Tk*(3.0495032E-7-Tk *3 . 4999996E-I0 ) ) 

2400 A2-< Nu5-2*Nu2+Nut )/200 

2405 A I -< Nu2-NuI-940*A2 )/ 10 

2410 A0-Nul-42*AI.-I764*A2 

2415 Nu-A0+Egr ♦ < A1 *Egr *A2 ) 

2420 RETURN Nu 

2425 FNENO 

2430 OEF FNCpeg( Tc , Egr ) 

2435 i RANGE OF VALIDITY 0 TO 20 OEG C 

2440 Tk=T c +273.15 

2445 Cpl-1 .670i5S0E+3+Tk*6.3 

2450 Cp2-1 .47481 25E+3+-Tk *6.25 

24S5 Cp5-9 . 5800S00E+2 + Tk *7 . 3 

2450 A2-< Cp3-2*Cp2+Cp l )/ 200 

24GS A1 -< Cp2-Cpl-900*A2 )/I0 

2470 A0«Cpl-40*AL-L600*A2 

2475 Cp-A0+E<3r* < A 1 +Egr * A2 ) 

2480 RETURN Cp 

2485 FNENO 

2490 OEF FNRhoeg< T , Egr > 

2495 Rol-1 .0607093E+3-T*< 3 . 703 1233E- 1 + T *4.0837 183E-3 ) 
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2500 Ro2-1.0748272E+3-T*<4.4266l95E-i+T*4„0939706E-3> 

2505 Ro3 a 1 . 0885934E+3-T ♦ ( 5. 73556 53E- 1 +T«6 . 1291 405E-3 ) 

2510 A2-(Ro3-2*Ro2*Rol )/200 

2515 Al-(Ro2-Rol-900*A2 )/ 10 

2520 A0-Rol-40*A1-1600*A2 

2S2S Ro-A0+Egr« < A1 +Egr*A2 ) 

2530 RETURN Ro 

2535 FNENO 

2S40 OEF FNPreg( T , Egr ) 

2545 Pr »FNCp eg < T , Egr )«FNNueg( T , Egr )«FNRoeg( T , Egr )/FNKeg( T f Egr ) 

25S0 RETURN Pr 

2555 FNEND 

2SS0 OEF FNKeg< Tc , Egr ) 

25S5 i RANGE OF VALIDITY : -20 TO 20 OEG C 

2570 Tk-Tc+273.15 

2575 Kl-2.2824708E-l+Tk •( S.S989236E-4+Tk«3 .5714286E-7 ) 

2530 K2*2.5846616E-l+TM(2.397357lE-4+Tk«7. 1423571E-7 ) 

2S85 K3-3.2133932E-1-Tk*< 3‘.0042357E-4-Tk* 1 . 4285714F.-6 ) ' 

2590 A2»< K3-2*K2+K 1 )/200 

2S9S Al-( K2-K l -900*A2 )/ 10 

2600 A0=K l -40^A 1- 1600 # A2 

2605 K-A0+Egr*( Al+Egr*A2 ) 

2610 RETURN K 

2615 FNENO 

2620 OEF FNTanh< X ) 

2625 P-EXP(X) 

2630 Q* 1/P 

2635 Tanh-< P-Q )/( P+Q ) 

2640 RETURN Tanh 

2645 FNENO 

2650 OEF FNTvav(U) 

2655 COM /Cc/ C( 7 ) 

2660 T-C<0) 

2665 FOR I-i TO 7 

2670 T-T+C( I >*V*I 

2675 NEXT I 

2630 RETURN T 

2635 FNENO- 

2690 OEF FN8eta(T) 

2695 Rop-FNRho<T+ . 1 ) 

2700 Ron»FNRho< T- . 1 ) 

2705 8eta--2/(Rop + RorO*< Rop-Ron ) / . 2 

2710 RETURN 8e ta 

2715 FNENO 

2720 OEF FNPsat ( Tc ) 

2725 » 0 TO 30 deg F CURVE FIT OF Psat 

2730 Tf *1 . 8 # Tc+32 

2735 Pa=»5 . 94S525 + T f • < . 1 5352032+T f * < 1 . 4340963E-3+T f *9 . 6 15067 IE-6 ) ) 

2740 Pg-Pa-14.7 

2745 IF Pg/0 THEN i + -PSIG,— in Hg 

2750 Paat-Pg 

2755 ELSE 

2760 Paat*Pg«29.92/14.7 

2765 ENO IF 

2770 RETURN Paat 

277 % 5 FNENO 

2905 SU8 Purge 

2910 SEEP 

2915 INPUT “ENTER FILE NAME TO 9E OELETED** ,F i le$ 

2920 PURGE Files 
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2925 

2930 

2935 

2940 

2945 

2950 

2955 

2960 

296S 

2970 

2975 

2976 

2977 

2980 

2985 

2990 

2995 

3000 

300S 

3010 

3015 

3020 

3025 

3030 

3035 

3040 

304S 

3046 

3050 

30S1 

3052 

3053 

3054 

30S5 

3056 

3057 

3058 

30S9 

3060 

3062 

3063 

3065 

3070 

3075 

3080 

3085 

3090 

3095 

3100 

3105 

3110 

3115 

3120 

3125 

3130 

3135 

3140 

3U5 

3150 

3155 



SOTO 2910 
SU8END 
SUB Wilson 



COM /Will/ Daa< 4 ) ,Dia< 4 ) ,Kna< 4 ) , lac t ,Qroot ( 4 ) , Ld< 3 ) 
COM /Wil2/ De 1 ta , Isa t , Nset s , Hod , C ia( 4 ) , A Ipaa < 3 ) 

OIM Enf< 17 ) ,T( 17 ) ,Xa< 100 ) , Ya( 100 ) ,Pc< 4 ) ,Tp< 4 ) 

8EEP 

INPUT "PLEASE RE-ENTER NAME OF FILE - , File$ 

ASSIGN SFiie TO File* 

INPUT ‘0-SMOOTH COPPER TU8E, l -CUN I FIN, 2-K0R0DENSE, 
8EEP 



INPUT "0-NO INSERT, 1- TWISTED TAPE INSERT, 2-HITRAN1 
BEEP 

INPUT "GIVE A NAME FOR XY FILE",Xy$ 

CREATE BOAT Xy$,20 
ASSIGN 9Xy TO XyS 
L * 1 .21 92 



Oo-Ooa( I tube ) 



Di-Qia< I tube ) 



K*-Kma< Itube ) 

Ax-PI*0i A 2/4 ! Cross-sectional area 



Ao-PI*Do*L 

Rn-Oo*LOG< Oo/Di >/<2*Km) 



3-FIN TITAN* , Itube 
3-HITRAN2", Insert 



Initial values 
Tf-Tsat 
Alpa- .65S 
IF Insert-0 THEN 
Ci-1 .0 • 

ENO IF 

IF Insert-1 THEN 

Ci-5.172 

ENO IF 

IF Insert-2 THEN 
Cl-. 226 
ENO IF 

IF Insert-3 THEN 

Ci-. 063 

ENO IF 

G-9.31 

Ibeg-0 

Iend-0 'CHANGE TO 4 , IF FIUE TU8ES IN BUNOLE 
IF Iact ' 10 THEN 
Ibeg- lac t 
lend- lac t 
END IF 

FOR I-Ibeg TO lend 
S/s -0 
Sy -0 
S a 5 — 0 
5xy-0 
Jset -0 

ASSIGN ©File TO F l 1 eS 

ENTER 9File.0tg$, Itube,Egrat ,0dl ,0d2 , 0d3 , Dd4 , 0d5 
ENTER SFile;Pc< * ) ,Enf < ♦ ), Tp< • ) 

FOR J-0 TO 17 

T< J )»FNT v sv< Ewf < J ) ) 

NEXT J 

Tvap »< T< 0 >+T < 1 >+T< 2 ) )/3 
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3160 T l iq-( T( 3 ) + T( 4 ) )/2 

3151 Tiiql-T(S) 

3165 IF Isat-0 THEN 

3170 Tsat-Tliql 

3175 ELSE 

3180 Tsat -Tvap 

3185 ENO IF 

3190 ' Grad-FNGrad(T< 1+6 ) ) 

3191 Tout »< T( 2*1 + 10 >+T( 2 + 1 + 11 ) )/2 

3195 Del t-Tout-T < 1+6 ) 

3200 Tavg-T< 1+6 )+0elt*.5 

3201 ! IF I-Iend THEN 

3202 ! Oeit-T< 10)-T<9> 

3204 ! Tavg*( T(10)-T<9)>/2 

3205 ! ENO IF 

3206 I 

3210 ! Water /Ethy lane Glycol Mixture Properties 

3215 Rhoeg-FNRhoeg( Tavg t Egrat ) 

3220 Nueg-FNNueg( Tavg , Egrat ) 

3225 Mueg-Nueg+Rhoeg 

3230 Cpeg-FNCpeg( Tavg f Egrat ) 

3235 Keg-FNKeg( Tavg f Egr3 t ) 

3240 Preg-Cpeg+Mueg/Keg 

3245 ! 

3250 Mdot -FNFncal ( I , T( I +6 ),Pc< I )) 

3255 Veg-Mdot / ( Rhoeg* Ax ) 

3260 Reeg-Veg+Oi/Nueg 

3265 Res*4+Mdot/< Mueg+( PI +0i-4+0el ta ) ) 

3270 Qdot-Mdot *Cpeg+0ei t 

3275 Qdp-Qdot/Ao 

3280 Lntd-Qeit/LQG< <Tsat-T< 1+6 ) )/( Tsat-T< 1+6 )-0elt ) ) 

3285 Uo-Qdp/Lntd 

3286 IF Insert-0 THEN 

3287 Bbi- . 0668*( Oi/L ) + Reeg*Preg 

3298 Bb2-l+.04*((Qi/L )*Reeg*Preg ) A . 6666 

3289 Onega-3 . 66 + ( 8b 1 /Bb2 ) 

3290 ENO IF 

3291 IF Insert-1 THEN 

3293 Omega-< 1 +5 . 484E-3+Preg* . 7 + < Res/Hod ) * 1 . 25 ) ~ .5 

3294 ENO IF 

3295 IF Insert-2 THEN 

3296 Onega- < Reeg* .65 ) *( Pr eg'* .46 ) 

3297 ENO IF 

3298 IF Insert-3 THEN 

3299 Qnega-( Reeg* . 76 ) * ( Preg* . 46 ) 

3300 ENO IF 

3302 i 

3303 1 R-U4/R1 13 Properties 

3305 Hf g*FNHf g( Tsat ) 

3310 Kf-FNK(Tf) 

3315 Rhof-FNRho( T f ) 

3320 Muf-FNMu(Tf) 

3325 ' 

3330 F-<Kr3 + Rhor2+G*Hfg/(Muf*0o + Qdp ) r . 33333 

3335 Ho-Alpa+F 

3340 Two-Tsa t-Qdp/Ho 

3345 Tf-rsat/3+2»Tuo/3 

3350 Y-( 1/Uo-Rn )*F 

3355 X=Qo*F/ ( Keg*0mega ) 

3360 i PRINT " OMEGA- " ; Onega ; " F» , F . * X- “ . X ; Y - " ; Y 
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3365 

3370 

3375 

3380 

3385 

3390 

3395 

3400 

3405 

3410 

3415 

3420 

3425 

3430 

3435 

3440 

3445 

3450 

3455 

3460 

3465 

3470 

3475 

3480 

3485 

3490 

3495 

3500 

3505 

3510 

3515 

3520 

3525 

3530 

3535 

8270 

8280 

8290 

8300 

8310 

8311 

3320 

8470 

8471 

8480 

8490 

3500 

8510 

8520 

8530 

8540 

3550 



i INEFFICIENT (MOOIFY LATER) 



Xa< Jaet )*X 
Ya< J set >-f 
Sx *Sx + X 
Sy-Sy+Y 
Sx s*Sx S + X *X 
Sxy*Sxy+X* Y 
Jset«Jset*L 
IF Jset<Nsets THEN 3135 
ASSIGN @Fxle TO * 

Slope*(Naets*Sxy-Sx*Sy )/<Naets # Sx3-Sx*2 ) 

Intcpt*< Sy-Slope*Sx )/Nsets 

Cic-1/Siope 

Alpac-1/ Intcpt 

Cerr-ABS< (Ci-Cic )/Cxc ) 

Aerr*A8S< < Alpac-Alpa )/Alpac ) 

IF Cerr> . 00 1 OR Aerr^.001 THEN 
Alpa*< Alpa+Alpac >♦ .5 
Cx-< Cx+Cxc )* .5 

PRINT *CIC*“ ; Cxc ; “ALFA*" ; Alpa 
GOTO 3100 
ENO IF 
8EEP 
BEEP 

PRINTER IS 1 

PRINT -CIC- - ;Cxc; “ALPA-- ; Alpa • 

Cxa< Insert )*Cxc 
Alpaa< I )*Alpac 
PRINTER IS 701 
FOR J*0 TO Nsets-1 

OUTPUT 9Xy,Xa< J),Ya< J) 

NEXT J 

PRINTER IS 1 
NEXT I 

ASSIGN 9Xy TO • 

SUBENO 

OEF FNFrcca 1 ( T ,Pc ) 

Mdot 15* . 0019*Pc- . 0098 
' fldo t0* . 00206* Pc - . 004 i 4 
Mdot 24* . 00225 ♦ Pc- . 00329 
IF T<0 THEN 

Mdt *< ( T+15 )*< Mdot0-Mdot 15 >/ 15 >+Mdot IS 
ELSE 

Mdt *< T *< Mdot24-Mdot0 )/24 )+Mdot0 
ENO IF 
RETURN Mdt 
FNENO 

OEF FNCosh( X ) 

P*EXP< X ) 

Q*EXP( -X ) 

Cosh* . 5* < P+Q ) 

RETURN Cosh 
FNENO 
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